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ABSTRAGT 

Integrating our knowledge of star formation traced by observations at different wavelengths is es¬ 
sential for correctly interpreting and comparing star formation activity in a variety of systems and 
environments. This study compares extinction corrected integrated ultraviolet (UV) emission from re¬ 
solved galaxies with color-magnitude diagram (CMD) based star formation rates (SFRs) derived from 
resolved stellar populations and CMD fitting techniques in 19 nearby starburst and post-starburst 
dwarf galaxies. The datasets are from the panchromatic STARBurst IRregular Dwarf Survey (STAR- 
BIRDS) and include deep legacy GALEX UV imaging, HST optical imaging, and Spitzer MIPS 
imaging. For the majority of the sample, the integrated near UV fluxes predicted from the GMD- 
based SFRs - using four different models - agree with the measured, extinction corrected, integrated 
near UV fluxes from GALEX images, but the far UV predicted fluxes do not. Further, we find a 
systematic deviation between the SFRs based on integrated far UV luminosities and existing scaling 
relations, and the SFRs based on the resolved stellar populations. This offset is not driven by different 
star formation timescales, variations in SFRs, UV attenuation, nor stochastic effects. This first com¬ 
parison between GMD-based SFRs and an integrated FUV emission SFR indicator suggests that the 
most likely cause of the discrepancy is the theoretical FUV-SFR calibration from stellar evolutionary 
libraries and/or stellar atmospheric models. We present an empirical calibration of the FUV-based 
SFR relation for dwarf galaxies, with uncertainties, which is ^53% larger than previous relations. 
Subject headings: galaxies: starburst - galaxies: dwarf - galaxies: star formation - galaxies: evolution) 


1. STAR FORMATION RATE MEASUREMENTS 

Measurements of the star formation histories (SFHs) 
of galaxies are a key to understanding the evolution 
of galaxies as a function of redshift. Qualitatively, the 
SFH can be deduced by classifying a galaxy as hosting 
an older, red population, a younger blue population, or 
both. Quantitatively, the SFH often reduces to a sin¬ 
gle measurable star formation rate (SFR) occurring in 
a galaxy at a particular redshift. Significant effort has 
been put into calculating such SFRs using a variety of 
techniques. One of the simplest methods is based on 
measuring the integrated light of a galaxy in a given 
wavelength regime (i.e., Ha, UV, or infrared (IR) wave¬ 
lengths), and converting this flux to a SFR using theoret¬ 
ically derived sc a hng relations (e.g.. iMadau et al.lllQ^ 
Keniaicuttl Il998t iGalzetti et al.l I2007L 120101: iHao et al l 
201lHMurphv et al.ll2011l) . These SFR scaling relations 
are based on assumptions about the origin of the in¬ 
tegrated light. For example, UV emission is assumed 
to originate from stellar populations with a well-defined 
range in age (and stellar mass), Ha emission is assumed 
to originate from gas ionized by massive stars younger 
than ^ 10 Myr, and light in the mid-IR regime is as¬ 
sumed to be emitted by dust that has been primarily 
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heated by relatively young stars. 

Galculating SFRs using simple scaling relations is most 
robust in the UV regime. In contrast to the Ha and 
IR emission that originates from gas and dust heated 
by star formation, the UV emission originates from the 
photospheres of young stars of intermediate and high 
mass (M> 3 Mq) and, thus, is a direct tracer of the 
recent star formation. SFRs derived from integrated 
emission measurements are dependent on the star for¬ 
mation timescales probed by a particular wavelength 
regime. The stars responsible for producing the major¬ 
ity of the UV flux have a lifespan of ~100 Myr, ver¬ 
sus ^5 — 10 Myr for the more massive stars (M> 17 
Mq) responsible for ionizing hydrogen and producing the 
secondary Ha emission. Galculation of SFRs based on 
integrated UV fluxes have been shown to be more ac¬ 
curate at lower values than SFRs based on integrated 
Ha fluxes, which can be suppressed due to the incom¬ 
plete sampling of the high mass end of the initial mass 
function (IMF) at low SFRs (e.g.. iBoselli et al.l 120091 
iGoddard et al.l [2?)T(il : iKoda et al.l 12012( 1. Similarly, tem¬ 
poral stochasticity can impact the interpretation of inte¬ 
grated UV e mission measur e ment s , although to a les ser 
extent (e.g.. iLee et al.l[200^ l201ll : IWeisz et al.ll201^ as 
the far ultraviolet (FUV) emission is produced by a wider 
range of stellar masses than the Ha emission. 

However, interpretation of the FUV emission can be 
complicated. Recent models of star formation (SF) in 
low-mass galaxies have found that different factors such 
as the total mass of a galaxy, the degree of clustering, and 
the overall level of activity may i mpact the total FUV 
emission produced at a giv en SFR (|Fumagalli et al.l[2011l : 
Ida Silva et al.l[20T^l20i^ . Further, the theoretically de¬ 
fined scaling relations are based on the assumption that 
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the UV emission originates in a system with a constant 
SFR. While this may be a valid approximation for spiral 
galaxies with a large number of independent star form¬ 
ing regions, this assumption can break down on smaller 
spatial scales or in low-mass systems. Similar to Ha emis¬ 
sion, UV emission can be significantly affected by extinc¬ 
tion, which can result in large uncertainties in measured 
SFRs. Best practices to accurately account for this at¬ 
tenuation require measurements of the energy absorbed 
by the dust and re-emitted at infrared wavelengths. Em¬ 
pirical studies using IR emission to trace dust content 
and approximate energy balancing schemes to convert 
this emission to an equivalent UV luminosi ty have im¬ 
proved corrections for UV attenuation (e.g.. feuat et al.l 
l2005t iCalzetti et al.l[200^ iKennicutt et al.ll2009ll . 

The Galaxy Evol utionary Explorer (GALEX) mission 
(jMartin et al.l[2?)?)^ surveyed two-thirds of the sky in the 
UV, greatly benefiting extragalactic studies in the UV. 
GALEX observations have been adept at characterizing 
star formation in low density gas environments, such as 
the outer disks of spiral and dwarf Irregular galaxies 
(|Thilker et al.ll200^ . providing a new means for tracing 
recent star formation in more quiescent environments. 
The ease of measuring the total integrated UV emission 
from GALEX images for a significant number of galaxies 
has made accurate calibration of UV-based SFR scaling 
relations even more important. 

Theoretical SFR relations based on scaling integrated 
UV emission measurements are derived and calibrated 
from synthetic spectral models of single- a ged stellar pop¬ 
ulatio ns (e.g., iBruzual fc Gharlo^I1993I iLeitherer et al.l 
B999I1 . They are based on the premise that if a SFR has 
been constant over a ~ 100 Myr timescale, the birth of 
stars that significantly contribute to the FUV emission is 
in equilibrium with the death of similar UV contributors. 
Thus, the UV emission is proportional to the SFR itself. 
UV scaling relations are widely used because of their rel¬ 
ative simplicity and because they ha ve been consistently 
calib rated against each other (e.g., IKennicutt fc Evansl 
I 2 OI 2 L and references therein). The most widely used UV 
scaling relation, presented in iKennic uttI (Il9 98fl. is based 
on the theoretical results from iMadan et al.l (jl998[l over 
the wavelength range 1500-2800 A. Adjustments to the 
scaling relation have been made using the stellar evo¬ 
lution models of the Starburst99 code (jLeitherer et al.l 
I 1999 H . In one study, the scaling fa ctor was raised by 14% 
for th e same IMF and metallicity (jIglesias-Paramo et al.l 
I 2 OO 6 II . but it is unclear what caused the increase. In 
a separate study, the scaling factor was lowered by 
10% bas ed on the lower met allicity values of the galaxy 
sample (iHunter et al.l 12010(1 . More recently, the the¬ 
oretical relation was re-calibrated using up dated stel¬ 
lar e v olution libraries and spectral models (jHao et al.l 
I 2 OIII: [Murphy etldl 20Ilh with the Starburst99 code 
(|Leitherer et al.lll999l : iVazauez fc Leithereil 1200511 . The 
resul ting revision is 5% lower than the original value 
from iKeimicutd (|I998ll for the same IMF. 

The FUV-based SFR scaling relation was empirically 
calib rated to be 30% l ower than the theoretical relation 
from IKennicuttI (jI998[l for the GALEX FUV bandpass 
using spectral energy distribution fitting of integrated 
light measurements for a large sample of high-mass galax - 
ies with a mean metallicity of 0.8 Zq (|Salim et al.ll2007ll . 


This study reported that calibration of this GALEX 
FUV-based SFR scaling relation was difficult for lower 
mass systems as proper treatment of the UV attenuation 
using only FUV and near ultraviolet (NUV) magnitudes 
was unreliable at lower SFRs. Further, the empirical 
calibration was derived from a sample of galaxies with 
a variety of SFHs. Thus, the assumption of constant 
star formation over the lifetime of UV emitting stars is 
compromised. 

Despite both the theoretical and empirical adjust¬ 
ments, what is still lacking is an independent, empiri¬ 
cal calibration of the UV scaling relation in a sample 
with well-measured and well-constrained star formation 
properties outside of the synthesis model approach. A 
different approach to measuring SFRs in nearby galaxies 
is to derive star formation histories (SFHs , or SFRs(t,Z)) 
from resolved st e llar populations (e.g., iDo luhinl 119971 
Tosi et aP 119891: iTolstov fc Sahal 119961: iGallart et al.l 


19961 : iHoltzman et al.l 11991^ iHarris k, Zaritskvl 1200111 . 

This method uses GMD fitting techniques to match ob¬ 
served CMDs to composite synthetic stellar populations 
generated from stellar evolution libraries. The best-fit 
modeled GMD represents the most likely SFH of the 
system. While deriving SFRs at the earliest epochs us¬ 
ing this technique requires deep optical imaging reaching 
below the oldest main sequence turn-off, recent SFRs 
can be robustly derived from well-populated CMDs of 
less photometric depth which include significant num¬ 
bers of upper main sequence and blue hel ium burning 
stars (lDofohinll2002l: iMcQuinn et al.l[2010bll . 

SFHs derived from resolved stellar populations pro¬ 
vide an independent measure of the SF activity which 
can be compared with SFRs from integrated UV mea¬ 
surements. In low-mass systems, which may be vulner¬ 
able to stochastic effects, the resolved stars also provide 
constraints on such effects. Such a comparison requires 
a sample of nearby galaxies with both significant lev¬ 
els of recent star formation and low levels of extinc¬ 
tion. Thus, nearby starburst dwarf galaxies make an 
ideal sample. By definition, starbursts have elevated 
levels of recent star formation activity, making them 
UV bright systems. Previ ously thought to have dura¬ 
tions of only ^ 10 Myr (e.g. . iSchaerer. Contini. fc KunthI 
19991: iMas-Hesse fc KunthI Il999t Thornlev et al. ^OOl : 


Tremonti et al.l l2001h . IMcQuinn et al.l (|2009L I2010M 
have argued that starbursts in dwarfs are a longer-lived 
phenomenon, lasting > 100 Myr, based on a comparison 
of their recent star-formation activity with their histor¬ 
ical averages derived from resolved stellar populations 
GMD-fitting techniques. From a companion study of the 
same data, the spatial extent of the star-formation ac¬ 
tivity was observed to range from centrally concentrated 
to widely distributed (IMcQuinn et al.lF201^ . The longer 
lived, often spatially distributed starbursts approximate 
a constant SFR over UV timescales - an assumption that 
underlies the SFR scaling relations. Dwarf galaxies also 
have lower metallicities than more massive systems, with 
correspondingly lower UV attenuation from dust. This 
reduces the uncertainties and complexities in applying 
extinction corrections to the UV emission measurements. 

In this study, we have undertaken a comprehensive 
comparison of the integrated UV emission from star for¬ 
mation based on deep GALEX imaging with the re¬ 
solved stellar populations imaged with the Hubble Space 
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Telescope (HST) in a sample of nearby starburst and 
post-starburst dwarf galaxies. These observations are 
part of the panchromatic STARBurst IRregular Dwarf 
Survey (STARBIRDS ) presented in a companion paper 
(jMcQuinn et al.ll2015tl . and have been used in previous 
studies to me asure the characteristics of th e star burst 
phenomenon (|McQuinn et'ani2009ll2010allbl . I2012tl . Al¬ 
though internal extinction was expected to be low at the 
metallicities of the sample, we applied the current best 
practice to correct for UV extinction using the infrared 
flux measured from Spitzer Space Telescope MIPS imag¬ 
ing. The data are discussed in Section 2. In Sections 3 
and 4 of this study, we test our understanding of the 
connection between optical and UV emission from stellar 
populations. Using four different models and the SFHs 
reconstructed from the resolved stellar populations, we 
predict the NUV and FUV fluxes from the galaxies, and 
compare these values with the extinction corrected, ob¬ 
served values of UV emission. This comparison provides 
a critical test for our understanding of the optical and 
UV emission in galaxies because it uses an independent 
measurement of the star formation from the resolved stel¬ 
lar populations to model the UV flux and then compares 
this modeled flux with the measured emission from ob¬ 
servations. Moreover, it provides a check on a number 
of the parameters including whether or not the SFHs are 
accurately describing the star formation activity in the 
galaxies, whether stochastic effects are dominating the 
results in low-mass galaxies, and whether the UV emis¬ 
sion has been accurately corrected for extinction. In Sec¬ 
tion 5, we compare the SFHs over timescales applicable 
to UV bright stars from resolved stellar populations with 
SFRs based on integrated FUV fluxes. In Section 6 we 
present an empirical calibration of the integrated FUV 
emission SFR scaling relation for dwarf galaxies based 
on the resolved stellar populations and in Section 7 we 
discuss possible sources of the dispersion in SFRs based 
on integrated FUV emission. Our conclusions are sum¬ 
marized in Section 8. In the Appendix, we provide a 
detailed discussion of the possible role observational un¬ 
certainties have on a comparison between predicted and 
measured UV fluxes. 

2. THE GALAXY SAMPLE AND OBSERVATIONS 

Table [T] lists the 19 nearby (D< 6 Mpc) starburst and 
post-starburst dwarf galaxies in the STARBIRDS sam¬ 
ple and their various physical properties. The sample 
was adopted from previous studies deriving the SFHs of 
nearby starburst dwarf galaxies using HST F606W (or 
F555W ) and F814W archival imaging (jMcOuinn et al.l 
l2010ajl . The sample was defined based on data from the 
HST archive that met the following criteria: First, the 
galaxies had to be close enough that the HST imaging in¬ 
struments could resolve the stellar populations. Second, 
both V-band (i.e., F555W or F606W) and I-band (i.e., 
F814W) observations were available for each system. 
Third, the I-band photometric depth reached a minimum 
of ~ 2 mag below the tip of the red giant branch. The 
galaxies span a wide parameter space in SFR, metal- 
licity, and physical size, including both stronger, well- 
known starburst dwarf galaxies in the nearby universe 
(e.g., NGC 1569, NGG 4449, NGG 5253), smaller, less 
well-known starbursts (e.g., NGG 4068, UGG 6456), and 
post-starburst galaxies (Antlia, UGG 9128, NGG 4163, 



Figure 1. Representative of our sample, GALEX NUV and FUV 
images of NGC 4068 are shown in the left and right panels re¬ 
spectively with North up and East left. The field of view matches 
the footprint of the HST ACS images, and excludes the chip gap, 
which facilitates a direct comparison of the star formation imaged 
by the different instruments. Typical of the sample, the image of 
NGC 4068 shows several knots of UV emission tracing the star 
forming regions, as well as fairly ubiquitous diffuse UV emission. 
The NUV image is significantly deeper than the FUV image and 
traces more of the lower surface brightness, diffuse UV emission. 

NGG 6789, NGG 625). Note, however, that it is not a 
volume limited sample. A more detailed discussion of 
the sample sele c tion a nd its limitations can be found in 
iMcQuinn et al.l (120151 see Section 2.1). 

The data consist of observations from the GALEX, 
HST, and Spitzer space telescopes. All reduced images 
are publically available through the STARBIRDS data 
archive available as high-level science products through 
the Mikulski Archive for Space Telescopes (MAST). The 
data reduction pr ocedures are describe d in detail in a 
companion paper (jMcQuinn et al.l[2015ll . Additional in¬ 
formation o n the processing o f the H ST observations can 
be found in IMcQuinn et aH (|2ni0afl . Below, we briefly 
summarize the data reduction process for each dataset. 

2 .1. GALEX UV Observations 

The GALEX instrument, with both FUV (A 1350 — 
1750 A) and NUV (A 1750 — 2750 A) b andpasses, is de¬ 
scribed in detail in iMartin et al.l (120051). and its i n -orbit 
perfo rmance and calibration by iMorrissev et al.l (|2005l 
l2007f ). The GALEX data used here include NUV and 
EUV images from multiple observing programs. The 
data were processed through the GALEX pipeline (v7.1) 
and all observations were co-added into one EUV and 
one NUV image per target. Figure [T] presents an example 
of the GALEX NUV and FUV images for the starburst 
galaxy NGG 4068, cropped to the field of view of HST 
observations of this system. In both images, UV bright 
knots from star forming regions are apparent, as well as 
extended, diffuse emission. The NUV image is consider¬ 
ably deeper than the FUV image and, thus, is sensitive to 
lower surface brightness features in the diffuse emission. 

The observations used to create the deeper, co-added 
tiles can have slightly varying fields of view. Therefore, 
the GALEX pipeline provides high resolution relative re¬ 
sponse (rrhr) maps which measure the depth of the ob¬ 
servations per pixel. In addition, the GALEX pipeline 
converts the GALEX satellite telemetry data (i.e., pho¬ 
ton lists) into calibrated intensity maps in counts per 
second and performs background subtraction on the co¬ 
added images. We used these background subtracted 
intensity maps for the lower surface brightness galaxies, 
but found that these images over-subtracted the back¬ 
ground for higher surface brightness regions in 5 galaxies 
(NGG 625, Ho H, NGG 6822, NGG 4214, NGG 4449) and 
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Figure 2. Representative of our sample, this HST composite 
image of NGC 4068 was created using the F606W image (blue), 
F814W image (red), and an averaged F606W and F814W image 
(green) with North up and East left. As the HST optical imaging 
provides a direct census of the resolved stellar populations, the HST 
footprint defined the field of view for comparing the star formation 
imaged by the different instruments. Apparent in the image are 
several blue star forming regions which correspond to the knots of 
UV emission seen in Figure^ An underlying redder, older stellar 
population is seen extending to larger radii. 

for one lower surface brightness galaxy affected by cirrus 
contamination (Antlia). In these cases, we performed a 
custom background subtraction. 

As we are interested in directly comparing the UV 
emission with the optical HST images, we cropped the 
larger, circular 1.25° GALEX field of view to the smaller 
footprint of the HST images. Finally, we custom masked 
the foreground and background contamination in the 
cropped, background subtracted intensity maps using the 
IRAFQ task APPHOT task imedit. The integrated 
NUV and FUV fluxes were measured in these post- 
masked images using the IRAF apert ure photometry 
functi onality. The fluxes are reported in iMcQuinn et al.l 
(|2015f) . and repeated here in Table [5] for convenience. 

2 .2. HST Optical Observations 

The HST observations are comprised of Advanced 
Camera for Surveys (ACS; iFord et ^ I1998D and Wide 
Field Planetary Camera 2 fWFPC2: iHoltzman et aP 
Il995f) archival images obtained with the F814W filter and 
either the F555W or F606W filter. The HST data were 
processed and cleaned by the standard HST pipeline. 
Figure [5] presents an example of the HST images for 
NGC 4068, the same galaxy shown in Figure [TJ This 
HST composite image shows the resolved stellar popula¬ 
tions including a young, blue stellar population and an 
underlying older, red stellar population. 

Photometr y was perforrn ed on the ACS images using 
DOLPHOT ((DolDhinll2000fl and on the WFPC2 images 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 


NGC 4068 
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Figure 3. Representative of our sample, this CMD is for 
NGC 4068 with average photometric uncertainties per magnitude 
bin shown at the left. The contours show the density of points 
where they would otherwise saturate the plot. The main branches 
of stellar evolution are noted. One of the distinguishing features 
in the CMDs of starburst galaxies is the notable population of red 
and blue helium burning stars and upper main sequence stars; both 
are evidence of significant recent star formation activity. 

using HSTphot (lDolDhinll2000D . Artificial star tests were 
performed to measure completeness due to the finite pho¬ 
tometric depth and crowding as well as the noise char¬ 
acteristics of the data using the same software. The re¬ 
sulting photometry star lists were uniformly processed, 
ensuring reliable comparison between ga laxies. More de¬ 
tails on the photometry can found in IMcQuinn et al.l 
(I2ni0all . 

In Figure[3l we show a CMD created from the photom¬ 
etry of the HST images for NGC 4068, which provides a 
direct census of the stellar populations. Evident in this 
representative CMD are distinct evolutionary branches 
of the stellar populations including the main sequence 
(MS), blue and red helium burning stars (BHeB and 
RHeB), the red giant branch (RGB), and asymptotic 
branch stars (AGB). Average uncertainties per magni¬ 
tude bin, including uncertainties measured from com¬ 
pleteness tests, are also shown. 

2.3. Spitzer IR Observations 

The Spitzer observations of the entire sample are com¬ 
prised of ar chival mid-IR obs e rvations from the MIPS 
instrument (|Rieke et al.l 12004 lEngelbracht et akl l2007f) 
in the 24, 70, and 160 ^m bandpasses. Emission at these 
mid-IR wavelengths traces the dust content heated by the 
stellar population and can be extrapolated to estimate 
the dust attenuation of the UV emission. Most of the ob¬ 
servations were obtained as part of eith er the Spitzer In¬ 
frare d Nearby Galaxy Su rvey (SINGS) (|Kennicutt et al.l 
[20(H IRegan et al.l [200l or the Local Volume Legacy 
iLyLl (|l)ale et all I2009H projects. These data were 
cleaned, background subtracted, and made publically 
available as value-added data products of the SINGS and 
LVL projects. Detailed descriptions of the data process¬ 
ing steps are provided by these legacy programs. Four 
galaxies (IC 4662, NGC 1569, NGC 6789, UGC 456) were 
not observed in either of these surveys. For these sys¬ 
tems, we used the MIPS observations available from the 
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Figure 4. Representative of our sample, MIPS 24 /im, 70 /im, and 
160 /im images of NGC 4068 are shown from left to right with North 
up and Each left. The MIPS images are used to estimate the total 
IR flux from the galaxy and the corresponding UV attenuation. 
The field of view matches the footprint of the HST ACS image, 
excluding the chip gap, enabling a direct areal match of the dust 
emission with the UV emission measurements. The MIPS images 
show a few knots of 24/im emission and more extended emission at 
the 70^m and 160/im wavelengths. The generally low level of IR 
emission is typical of this low metallicity sample. 

Spitzer Heritage Archive (SHA). These archival obser¬ 
vations covered only 2 of the 3 MIPS bandpasses, and 
thus provide only lower limits of the dust attenuation of 
UV emission. We cleaned and mosaicked the SHA ob- 
servations using the Spitzer-provided MOPEX software 
(iMakovoz fc Khanl 120051 : iMakovoz et al.l l2005f) . Back¬ 
ground subtraction was performed using an average flux 
value from relatively blank regions of the observations. 

The MIPS images were treated in the same way as 
the GALEX images. First, the MIPS observations were 
cropped to match the field of view of the HST im¬ 
ages. Second, we custom masked the foreground and 
background contamination in the same manner as the 
GALEX UV images. Finally, the integrated mid-IR 
fluxes were measured in the post-masked images using 
the IRAF task APPH QT. The fluxes are reported in 
iMcQninn et al.l ()2015[1 . and repeated here in Table [2] for 
convenience. 

Figure H] presents an example of the MIPS images at 
24 /rm, 70 /rm, and 160 /rm for NGC 4068, the same 
galaxy shown in Figures [THHl Typical of most of the 
sample, the MIPS images show only a low level of emis¬ 
sion, with little emission detected at 160 /rm. Since the 
majority of the sample has relatively low values of metal¬ 
licity, the mid-infrared emission from dust heated by star 
formation is also relatively low. The MIPS data are used 
to estimate the total IR flux and corresponding UV at¬ 
tenuation. 

3. DERIVING SFHS FROM RESOLVED STARS, 
PREDICTING UV FLUXES, AND CALCULATING SFRS 
FROM INTEGRATED FUV FLUXES USING EXISTING 
SCALING RELATIONS 

Our calculations follow three stages. First, we describe 
the derivation of the SFHs from the resolved stellar popu¬ 
lations. Second, we use four models in conjunction with 
the derived SFHs, to predict the integrated NUV and 
FUV fluxes in our sample. Third, we convert the UV 
and MIPS fluxes to extinction corrected UV luminosi¬ 
ties and apply existing scaling relations to convert the 
FUV luminosities to SFRs. A discussion of these results, 
including comparisons between the measured and pre¬ 
dicted quantities, is presented in Sections 4 and 5. 

3.1. Reconstructing SFHs from Resolved Stellar 
Populations 


The SFHs were previo usly reconstructed by 
iMcQuinn et al.l (I2009L l2010all . Briefly, the SFR(t,Z) 
was derived b y applying a GMD fitting technique 
(jPolphinl l2002tl with stellar evolutionary models to the 
photometry and artificial star l ists. We used the Padova 
isochrones (jMarigo et al.l 1200^ because they cover the 
range in metallicity of our sample and include isochrones 
for the high mass stars that are of interest to us. The 
GMD fitting technique fits for distance, extinction, star 
formation history, and chemical enrichment history, with 
the metallicity constrained to increase with time (except 
for Antlia and NGG 2366 where the photometry reached 
a full magnitude be low the red clump). A Salpeter 
IMF (iSalpeterl II955I1 with mass limits of 0.1—120 Mq 
is assumed in the derivation. A single power-law IMF 
is sufficient for our purposes as the stars of interest are 
all high-mass, UV producing stars. The distances fit 
by the method were in good agreement with tip of the 
red giant branch (TRGB) distances from the literature, 
which provides a consistency check on the GMD fit. 

In Figure [5l we reprod uce the SFR(t) over th e ^ 0.5 
Gyr for NGG 4068 from IMcQuinn et al.l (j2010b[ l. Note 
the elevated SFR at recent times indicative of the star- 
burst and representative of the SFHs derived for the sam¬ 
ple. The uncertainties in the SFR(t) include both sys¬ 
tematic and statistical uncertainties derived from Monte 
Garlo simulations. The time binning reflects the inherent 
time resolution of the data. At older look-back times, the 
SFRs are not well constrained without photometry that 
reaches the oldest MS turn-off. However, at the more re¬ 
cent times shown here, the low photometric errors of the 
massive, young MS and HeB stars enable robust mea¬ 
surements of the SFRs with a m uch higher temp o ral res- 
olution. We refer the reader to IMcQuinn et al.l (|2010a[) 
for a full discussion of the SFHs, the time binning, and 
the calculation of un certainties. _ 

The SFRs(t) from IMcQuinn et al.l (|2010ar) were used 
to calculate the average SFRs over the past 100 Myr. 
These GMD-based SFRs are listed in Table p We chose 
this timescale because 90% of the FUV emission in galax¬ 
ies with a co nstant SFR origi nates from stars less than 
100 Myr old (jHao et al.l 1201 Ih . This temporal baseline 
changes to 200 Myr when considering NUV emission 
(iHao et al.ll2T)Illl . To conservatively estimate uncertain¬ 
ties, we report the standard deviation in the SFRs over 
the past 50, 100, and 150 Myr. In all cases, the stan¬ 
dard deviations are larger than the measured uncertain¬ 
ties from the SFH. The benefit of taking this approach 
is that it provides a measure of how constant the SFRs 
are in each galaxy over the timescale of interest. 

3.1.1. Special Cases 

As with any SFR measurement, a GMD-based SFR 
depends on the observations to accurately trace the star- 
formation activity. The high-resolving power of the 
HST and the application of artificial star tests means 
that GMD-fitting techniques can successfully measure 
the recent star-formation in galaxies out to several Mpc. 
However, we have identified four galaxies in our sam¬ 
ple which host young, massive clusters that are unre¬ 
solved with the HST imaging. Thus, the GMD-based 
SFRs may be under-reporting the true star-formation 
activity. Further, due to the significant number of very 
young stars (t< 10 Myr), the assumption of constant 
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Figure 5. The SFH of NGC 4068 over the last 600 Myr recon¬ 
structed from the r esolved stellar populations using the CMD fit- 
ting method MATCH IIDolDhinll200^ and stellar evolutionary models 
IIMarigo et al.l2n08ll . The time resolution improves at recent times, 
with time bins as short as 20 Myr. The SFRs are elevated over the 
past ~ 500 Myr, reflecting the starburst activity in this system. 
While the SFH of each galaxy is unique, SFRs elevated for more 
than 100 Myr is a common feature in the sample. 

star-formation over the UV timescale is compromised. 
Thus, the FUV-based SFRs may be over-reporting the 
true star-formation activity. 

Specifically, NGC 5253, NGC 4449, NGC 4214, and 
NGC 625 are known to host young (t< 10 Myr), mas¬ 
sive 10^ Mq) clusters. In NGC 5253, the dy¬ 
namical mass of th e central star cluster is es timated 
to be 3 X 10® Mq (|Rodriguez-R.ico et al.ll200^ . while 
numerous clusters with stellar mass of 10^ Mq have 
also been measured (|de Griis et al.l[20T^ . Similarly, in 
NGC 4449, the centr al super star clust er has a stellar 
mass of 2 X 10® Mq (|B5ker et al.l l200l[l with up to 20 
addit ional star clusters w ith masses between 10'’’ — 10® 
Mq (jAnnibali et al.ll2011fl. NGC 625 hosts a central¬ 
ized starburst (|Cannon et al.l|2003fl : the Ha luminosity 
of the central cluste r approaches that of 30 Doradus 
(|Skillman et al.l I2003D whose ste llar mass is estimated 
to be 4 X 10® M.;. (|Bosch et all [200^1 . NGC 4214 in¬ 
cludes multiple massive clusters in the cent ral region with 
mass es ranging from 2x 10^ — 5x 10® Mq (|Sollima et al.l 

Hall). 

Visual inspection of the HST optical images for these 
four galaxies used to derive the CMD-based SFRs con¬ 
firms knots of unresolved stars in these massive clus¬ 
ters. In unresolved clusters, stars that are not photo¬ 
metrically recovered because of blending do not factor 
into the CMD-based SFRs, but certainly contribute to 
the integrated UV emission measured in the GALEX 
images. Therefore, a comparison between CMD-based 
SFRs and FUV-SFRs will be biased. We include these 
outliers in the rest of the study for completeness, differ¬ 
entiating them from the rest of the sample with unique 
plot symbols. However, we do not include them in the 


calculation of subsequently reported values or in the SFR 
calibration. 

3.2. Predicting Integrated UV Fluxes from CMD-based 

SFRs 

The SFRs(t,Z) derived from the resolved stellar popu¬ 
lations can be used to model the present-day extinction- 
free UV luminosity expected from a galaxy, taking 
into account both the changes in SFR over time and 
the metallicity of each system. We used four mod¬ 
els with different approaches to predict the UV lu¬ 
minosities from the SFHs: (1) the synthetic spec¬ 
tral model PEGASE 2.(0 (|Fioc fc Rocca-Volmeran^ 
1199711 : (2) synthetic s tellar populations generated from 
MATCH (jDolphirJ 1200211 which stochastically populates a 
synthetic UV CMD based on an IMF and the input 
SFR and metallicity evolution; (3) the synthetic stel¬ 
lar populati on model Stochastical l y Lighting Up Galax - 
ies (SLUG: [Fumagalli et al.l [20111 Ida Silva et ani2012ll . 
which stochastically populates the IMF and includes ad¬ 
ditional temporal stochastic effects due to clustered star 
formation; a nd (4) the stellar populati on synthesis model 
GALAXEV (iBruzual fc CharlotllMia . 

3.2.1. Predicting UV Fluxes Using the Synthetic Spectral 
Model Pegase 

For the first method, we used the synthetic spectral 
model Pegase 2.0 to generate a full spectrum for each 
galaxy from our SFHs. The Pegase code builds a syn¬ 
thetic galaxy spectrum from a set of global parameters 
such as an initial gas mass, modeled stellar feedback pro¬ 
cesses, chemical enrichment history, and SFH. The model 
assumes a fully populated IMF and scales the spectrum 
by the overall mass of a galaxy. We selected a Salpeter 
IMF with the mass limits of 0.1 — 120 Mq, to match that 
used in deriving the SFHs. The use of a single power-law 
IMF is appropriate as the the stars of interest (both in 
generating the recent CMD-based SFRs and those gen¬ 
erating the UV flux) are higher mass stars not impacted 
by a double power-law IMF. The Pegase model includes 
stellar evolutionary tracks from the main sequence to the 
white dwarf stage. The older Padova stell ar evolution¬ 
ary isochrones are used (iBertelli et al.lll994ll without the 
updates from iMarigo et al.l (1200811 . Dust extinction is 
modeled using a radiative transfer code. The main input 
to the Pegase model is the SFR(t,Z) normalized by the 
total baryonic mass of the system. The baryonic mass 
normalization factor was estimated in two steps. First, 
we estimated the total stellar mass created in the field of 
view of the HST observations in each galaxy by integrat¬ 
ing the lifetime SFHs over all time intervals. Second, 
the total gas mass was estimated to be approximately 
equivalent to the stellar mass, as gas-rich d warf irregulars 
typic ally have gas mass fractions of 0.5 (iBegum et al.l 
l2008f) . Therefore, the normalization factor of total bary¬ 
onic mass was 2 x E (SFR(At)x At). 

In order to include an estimate of uncertainties, we 
generated synthetic spectra using 50 Monte Carlo simu¬ 
lations of the SFR(t,Z). The final predicted spectrum is 
the mean of these 50 spectra and the dispersion on the 

** PEGASE 2.0 is publically available via 
http: //www2.iap.fr/users/fioc/PEGASE.html 


























































Calibrating UV Star Formation Rates for Dwarf Galaxies from STARBIRDS 


7 



Figure 6. Synthetic spectrum of NGC 4068 generated using the 
Pegase model and the SFR{t,Z) derived from optical imaging. The 
uncertainties are based on Monte Carlo simulations of the SFHs. 
The NUV and FUV bandpasses are highlighted in grey. 

mean provides an estimation of the uncertainty on the 
spectral luminosity. The predicted integrated UV fluxes 
were calculated by summing the Pegase present day 
spectrum across the GALEX NUV and FUV bandpasses, 
taking into acc ount the GALEX sens itivity curves for the 
UV detectors (|Morrissev et al.l[2005ll . Figure [6] presents 
an example of the synthetic spectrum for NGC 4068 
with uncertainties based on the standard deviation from 
Monte Carlo simulations. Both the GALEX NUV and 
FUV bandpasses are highlighted in grey. 

3.2.2. Predicting UV Fluxes Using Synthetic Stellar 
Populations from MATCH 

For the second method, we generated synthetic stellar 
populations based on the SFRs(t,Z) using the fake rou¬ 
tine from the CMD fitti ng code MATCH and Padova stel¬ 
lar evolutionary models (|Marigo et al.l[2008l) . The stellar 
populations were generated assuming the distances listed 
in Table [U and a Salpeter IMF with the same mass lim¬ 
its used in deriving the SFHs. The IMF is randomly 
populated for the input SFR(t). This is an important 
stochastic effect to consider at low-SFRs as the upper- 
end of an IMF may not be fully populated, and, thus, a 
star forming region may be under-luminous for a given 
stellar mass (or SFR). Additionally, internal extinction 
was set to zero in the model. The output from MATCH is 
a list of synthetic magnitudes based on the Vega mag¬ 
nitude system with luminosities in the HST F555W or 
F606W (depending on the hlter used in the actual HST 
observations), and F814W filters, and the GALEX NUV 
and FUV filters, using the calibration of the isochrones 
for the GALEX bandpasses. The UV luminosities of the 
synthetic stellar populations were converted from Vega 
magnitudes to the AB magnitude system using the con- 
ve rsions for the GALE X bandpasses from the database 
in iMarigp et al.l (|2008ll : 

NUVab = NUVvega + 1-662 (1) 

FUVab = FUVvega + 2.128 (2) 

As a consistency check, synthetic CMDs generated in 
the HST filters can be compared with the observed 


Synthetic Optical/NIR CMD for NGC 4068 

-9 


-8 



F606W - F814W (mag) 


Figure 7. Synthetic CMD in the HST filters for NGC 4068 based 
on the SFH reconstructed from the optical imaging. The main dif¬ 
ference between this synthetic CMD and the observed CMD shown 
in Figure [ 3 ] is the absence of observational errors in the synthetic 
stellar population. While it is possible to include errors using arti¬ 
ficial star tests, this introduces additional uncertainties when pre¬ 
dicting the UV flux expected from the CMD. 


Synthetic UV CMD for NGC 4068 

-9- ■ 


-R 



FUV - NUV (mag) 


Figure 8. Synthetic CMD in the GALEX filters for NGG 4068 
based on the SFH reconstructed from the optical imaging. This 
CMD is the UV companion to the synthetic CMD in the HST 
filters shown in Figure [7| The majority of the UV emission is 
produced by the upper MS and blue HeB stars. 


GMDs. Figure [7] shows the synthetic CMD generated 
for NGC 4068 from the SFH, which can be compared 
with observed CMD shown in Figure [3] The main dif¬ 
ference between the two CMDs is due to the generation 
of the synthetic stellar populations without observational 
errors. Figure [5] shows the corresponding synthetic CMD 
for NGC 4068 in the GALEX filters. The predicted in¬ 
tegrated UV luminosities were determined from the syn¬ 
thetic stellar populations by summing the UV luminosi¬ 
ties in the GALEX filters. 
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3.2.3. Predicting UV Fluxes Using Synthetic Photometry 
from SLUG 

For the third method, we generated integrated UV 
fluxes based on the SFRs(t) at the present day metallicity 
using SLUG. SLUG is a modification of the Starburst99 
code which uses Monte Garlo realizations to account for 
the sampling stochasticity of the upper-end of the IMF, 
temporal stochasticity from the finite sampling of the 
SFH, and stochastic sampling of the initial cluster mass 
function (IGMF). The SLUG model creates and evolves 
a synthetic galaxy over a maximum of 1 Gyr and gen¬ 
erates the expected flux at different bandpasses, includ¬ 
ing the GALEX NUV and FUV bands. Fo r consistency, 
the P adova stellar evolutionary models of iMarigo et al.l 
(|2008f) were used and the same IMF slope was assumed 
but with mass limits between 0.1 and 100 M 0 . However, 
in contrast to the other models, the galaxies are mod¬ 
eled using values of their present day metallicity and 
no chemical evolution is included. This will not have 
a significant effect on our comparison as there is little 
chemical evolution expected over the last 1 Gyr. How¬ 
ever, the metallicity inputs for SLUG are approximate 
as the choice of isochrones is limited to 0.0004, 0.004, 
and 0.008 without any intermediary interpolation. This 
will have some effect on the predicted FUV fluxes as the 
FUV emission from stars varies with metallicity (e.g., 
iMadau fc Dickins^l2014H . We discuss the dependency 
of FUV emission on metallicity further in Section 7. 

The SLUG model was developed to probe the effects 
of stochastic variations that can occur in stellar popu¬ 
lations at low absolute SFRs. Similar to MATCH, SLUG 
randomly populates the upper-end of an IMF for a given 
input SFR. Additionally, SLUG can be fine-tuned to al¬ 
low for unclustered or clustered star formation. In this 
context, a ‘cluster’ is meant to include any type of stellar 
association, dynamically bound or not. For a clustered 
star formation scenario, SLUG randomly populates each 
individual cluster with the assumed IMF. Generally, this 
results in lower UV fluxes for clustered SF, particularly 
if the cluster masses are constrained to be low-mass clus¬ 
ters, as UV-bright, high-mass stars have a lower proba¬ 
bility of occurring in low-mass clusters. 

3.2.4. Predicting UV Fluxes Using Stellar Population 
Synthesis Models from GALAXEV 

For the fourth and final method, we used the stel¬ 
lar population synthesis model GALAXEV to generate 
a full spectrum for each galaxy from our CMD-based 
SFHs. Similar to Pegase, the GALAXEV code builds 
a synthetic spectrum of a galaxy based on the SFH but 
assumes a constant metallicity over the lifetime of the 
system. This will not have significant effect on our com¬ 
parison as there is little chemical evolution expected over 
the last 1 Gyr. However, similarly to the SLUG model, 
the metallicity inputs are approximate as the choice of 
isochrones is limited to 0.0004, 0.004, and 0.008 without 
any intermediary interpolation. For consistency with the 
other predictions, we assumed the same IMF slope. How¬ 
ever, the mass limits for a Salpeter IMF in GALAXEV 
is set from 0.1 and 100 M 0 (see below for a discussion 
on the impact of the IME mass limi ts). In addition, the 
Padova stellar evo lutionary models (IMarigo et al.ll2008D 
were used. Erom l^uzual &: Gharlotl ( 20031) . the older 


1994 Padova models are preferred over these because of 
differences in temperatures of the red giant branch stars. 
However, not only do these differences not impact the 
UV fluxes, but constancy over the choice of stellar evo¬ 
lution models in the predictions is preferred. No dust 
attenuation was included in the spectra generated using 
GALAXEV. 

3.2.5. Comparison of Inputs and Assumptions Between the 

Four Models 

Wherever possible, we used the same inputs and basic 
assumptions in all four models to predict the UV fluxes. 
As described above, all the models use the GMD-based 
SFRs(t) measured from the resolved stellar populations 
and the Padova stellar evolution libraries. Further, no 
reddening corrections were applied making the UV lu¬ 
minosities all intrinsic, dust-free luminosities. However, 
there are some differences between the models. The most 
important difference is in the choice of stellar atmosphere 
models used to transform each model’s results into a 
spectrum or flux measurement. We discuss these dif¬ 
ferences in detail below (Section 4). 

There are also smaller differences between in the mod¬ 
els. For example, MATCH applies a user-specified binary 
fraction of 35% to generate binary stars in the synthetic 
populations based on a fiat secondary distribution. In 
contrast, Pegase, SLUG, and GALAXEV do not con¬ 
sider the contribution of binaries to the overall galaxy 
flux or spectrum. A Salpeter IME with mass limits from 
0.1 to 120 M0 was assumed in the Pegase and MATCH 
models; SLUG and GALAXEV used a slightly smaller 
range of 0.1 to 100 M 0 . This will have a negligible ef¬ 
fect on our results as few, if any, of these higher mass 
stars are expected at these low SERs. We found that 
changing the IMF mass limits to 0.1—100 M 0 had less 
than a 5% increase in the predicted fluxes from the syn¬ 
thetic spectral model Pegase. Both Pegase and MATCH 
utilize the chemical enrichment history information de¬ 
rived from the stellar populations and interpolate be¬ 
tween isochrones, while SLUG and GALAXEV assume a 
constant metallicity over the past 1 Gyr approximated at 
one of three isochrones. As noted above, the enrichment 
history will likely have a negligible effect on our results 
as there is little expected chemical evolution in the last 1 
Gyr, but small differences can be introduced when using 
approximate metallicities. 

As mentioned above, the SLUG model can be altered 
from populating the IMF in an unclustered star forma¬ 
tion scenario to one with clustered star formation. We 
found the UV flux predictions from SLUG depended on 
whether the star formation was clustered or not, but were 
not as sensitive to the degree of clustering (i.e., there was 
an appreciable effect on the flux predictions between as¬ 
suming no stars formed in clusters versus 100 % of the 
stars formed in clusters, but there were only small differ¬ 
ences between assuming 50%, 70% or 100% of the stars 
formed in clusters). Additionally, the SLUG predictions 
depended on the range of assumed cluster masses. 

3.3. Calculating Extinction Corrected UV Luminosities 

The measured NUV and FUV fluxes in Table [5] were 
converted to luminosities in the AB magnitude sys¬ 
tem using the calibration and zeropoint magnitudes 
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(|Morrissev et al.l[2007li : 

NUV = -2.5 • logio (FIuxnuv 

/ 2.06 X 10"^® erg sec"^ cm-^ A”^) + 20.08 (3) 

FUV = -2.5 • logio (Fluxpuv 

/ 1.40 X 10"^® erg sec"^ cm'^ A”^) + 18.82 (4) 

These UV magnitudes are based solely on the UV flux 
measured in each galaxy and do not account for possi¬ 
ble UV extinction along the line of sight. Therefore, we 
estimated the attenuation expected at UV wavelengths. 
Although many of the galaxies in our sample are low 
metallicity dwarfs, and thus are expected to have low in¬ 
ternal dust content and correspondingly low values of UV 
extinction, eight systems have somewhat higher metallic- 
ities (i.e., 8.1 < log 12 -I- (0/H) < 8.4; see Table[T]). 

One way to estimate the amount of attenuation at UV 
wavelengths is to measure the total IR emission from 
the dust, which is primarily heated by UV emission from 
stars, and convert this IR flux to an equivalent UV lu¬ 
minosity. This method is an improvement over using 
the FUV-NUV color to estimate extinction, which re¬ 
lies on assumptions about the intrinsic color of the emit¬ 
ting stellar populations, dust extinction, g eometry of the 
system, and the effects o f scattering fe.g.. iGordon et aH 
120011 iMisselt et ^l2001fl . The relationship between at¬ 
tenuation of UV fluxes a nd mid-IR emission from dust 
was deriyed empirically by lBuat &: Xul (1199611 for a range 
of extinction mechanisms. More recently, the attenua¬ 
tion expected at the GALEX bandpasses extrapolated 
from a measurement of the total IR flux was deriyed 
using models of actiye star forming regions (|Buat et al.l 
1200511 . While this prescription cannot differentiate be¬ 
tween IR flux from dust heated by UV emitting stars 
and dust heated by stars emitting longwards of the UV, 
differences in results from more massive starburst galax¬ 
ies with larger populations of UV emitting stars and from 
normal star-forming galaxies are only of order 30% (e.g. , 


Buat et al 


Hao et al.l 


Buat et al 


ll99(^[Meurer et ^Il999l [Gordon et al.ll2000l : 
201111. We use the empirical relations from 


(|2005[1 derived from a sample of spiral and ir¬ 
regular galaxies that extend to an M^uv = 16 mag (AB 
system), the approximate lower limit of our sample. We 
reproduce the relations here in Equations [5] and El 

Anuv = -0.0495 • + 0.4718 • + 0.8998 • a; -H 0.2269, 


where x = log(FTiR/FNuv) 


( 5 ) 


Afuv = -0.0333 • -f 0.3522 ■ -f 1.1960 • y -f 0.4967, 

where y = log(FTjj^/Fpuy) (6) 

The total IR flux needed for the above equations can be 
estimated using the MIPS flux measurements follow ing 
the empirical prescription from iDale &: Heloul (|2002ll : 

Ftir = 1.559 • 1 / ■ A(24^m) -h 0.7686 • i/ ■ A(70^m) 

-f 1.347 ■j^-A(160/rm) (7) 

The MIPS fluxes in Table E] and the relationships in 
Equations m were used to correct the NUV and FUV 
luminosities for extinction. Table [3] lists the extinction 


estimates (Anuv (mag) and Apuv (mag)), and the ex¬ 
tinction corrected NUV and FUV luminosities. For com¬ 
pleteness, this correction was performed for the entire 
sample, regardless of metallicity. Two galaxies, Antlia 
and UGG 9128, are non-detections in all three MIPS im¬ 
ages; the UV luminosities were not adjusted for these sys¬ 
tems. Four galaxies (UGG 6456, IG 4662, NGG 1569, and 
NGG 6789) were not observed in all three MIPS chan¬ 
nels; their estimated extinction values are lower limits as 
noted in Table El 


3.4. Calculating SFRs from Existing Scaling Relations 

As discussed in Section 1, there are a number of the¬ 
oretical scaling relations to convert UV luminosities to 
SFRs. Both the earli er calibration of the FUV-based 
SFR scaling relation bv iKennicnttJ (1199811 and the r ecent 
updated calib r ation s from both iHao et al.l (|2011h and 
iMurphv et al.l (1201111 are based on the Starburst99 spec¬ 
tral m odels (|Leitherer et ^119991 1 Vazquez fc Leithererl 
1200511 . The Starburst99 code uses an isochrone synthesis 
technique to generate a synthetic spectrum for a galaxy 
based on a given set of user-specified parameters (for 
example, a constant SFR and constant metallicity over 
some time period). Briefly, Starburst99 uses a stellar 
evolutionary library and an assumed IMF to generate 
a Hertzprung-Russell diagram, and then applies stellar 
atmospheric mod els to points in the binned H ertzprung- 
Russell diagram (iLeitherer fc Heckmanlll99^ . The cu¬ 
mulative sum of the modeled spectra results in a syn¬ 
thetic galaxy spectrum. The monochromatic flux in the 
synthetic spectrum at an FUV wavelength (generally 
taken as 1500 A) is then used as a basis to calibrate the 
theoretical integrated FUV-based SFR scaling relation. 

We chos e to u se the following scaling relation from 
IHao et aP (120111 ) based on the updated models from 
Starburst99 to convert our extinction corrected, inte¬ 
grated UV luminosities into SFRs: 


SFR = 1.33 X 10-28 X L°, 


( 8 ) 


where SFR is in Mq yr ^ and L)) is the attenuation cor¬ 
rected FUV luminosity in ergs s”^ Hz”^. The FUV- 
based SFRs calculated in this way are listed in Table El 
The conversion factor in Equation El was established 
specifically for the GALEX FUV bandpass and assumes 
a Salpeter IMF with mass limits of 0.1 — 100 Mq. The 
upper mass limit of 100 Mq is lower than the 120 Mq 
limit used in deriving the SFHs and in two of the mod¬ 
els used for predicting the UV fluxes. We found through 
experimenting with different mass limits in the synthetic 
spectral model Pegase that this upper limit has < 5% 
effect on the results (see discu ssion above). The SFRs us¬ 
ing the scaling relation from IHao et al.l (|2011ll are ~5% 
lower than tho se calculat e d from the theoretical scal¬ 
ing relation by iKennicuttI (11998( 1 and 23% higher than 
the SFRs c a lculat ed using the empirical relation from 
I Salim et'ni ()2007[ 1. However, as these scaling relations 
differ by a linear offset, it is simple to consider the effect 
of this choice on our results (see Section 5). 


4. COMPARISON OF THE MEASURED AND PREDICTED 
UV FLUXES 

One of the critical tests to our understanding of star 
formation measurements is whether the measurements 

































































10 


McQuinn et al. 


across wavelength regimes agree. For this study, the dif¬ 
ferent assumptions used in determining the CMD-based 
and FUV-based SFRs and the calibration of the FUV 
scaling relation are embedded in the comparison. In this 
section, we simplify the comparison by eliminating the 
calibration of the FUV-based SFR. We compare only 
the UV fluxes modeled from the independently mea¬ 
sured CMD-based SFRs with the UV fluxes measured 
from GALEX images. This test therefore integrates the 
resolved stellar population measurements with the syn¬ 
thetic spectral models and stellar population synthesis 
models, without the assumptions in the FUV-SFR scal¬ 
ing relation. In the following section, we will expand 
our analysis and examine the SFR values from both the 
CMD analysis and the FUV scaling relations. 

4.1. NUV Ratios 

The ratios of the measured to predicted NUV fluxes 
from the four models are reported in Tableland shown 
as a function of B-band luminosities in Figure ID The 
FUV ratios are shown in Figure [TO] and are discussed sep¬ 
arately in Section 4.2. Each panel in Figure [9] presents 
the results from one of the four models as labelled. Star- 
burst galaxies are plotted as filled symbols and post- 
starburst galaxies are plotted as unfilled symbols. The 
thin solid lines represent a perfect match between the 
predicted and measured values, the dashed lines repre¬ 
sent the mean value assuming a log-normal distribution, 
and the shaded regions represent the standard deviation 
on the log-normal mean. For completeness, we include 
the four galaxies whose CMD-based SFRs were previ¬ 
ously identified as lower-limits (plotted as stars), but do 
not include their NUV flux ratios in the mean. As ex¬ 
pected, the predictions for these galaxies are biased low. 

From Figure El the NUV fluxes predicted by all four 
models agree with the measured fluxes for the majority of 
the sample. This is remarkable given the complex inputs 
inherent in this comparison. The mean ratios and stan¬ 
dard deviations based on a log-normal distribution range 
from X = 0.94,0- = 0.24 to a; = 1.22, o' = 0.16 as noted 
in each panel. The generally high dispersion is not un¬ 
expected. The modeled fluxes are based on the SFRs(t) 
measured from the resolved stellar populations. How¬ 
ever, even with the fine temporal resolution achievable 
from the STARBIRDS data on these galaxies at recent 
times, the SFRs(t) are not continuously sampled. Fluc¬ 
tuations in the SFR are known to occur i n the sample on 
timesc ales shorter than a few Myr (e.g., iMcQuinn et al.l 
l20I0bf) . This temporal stochasticity can occur more fre¬ 
quently at lower SFRs in low-mass galaxies as the SFRs 
are measured over a small number of individual star¬ 
forming events. The standard deviations between the 
predicted and measured flux values are consistent with 
the types of variations one would expect given the finite 
temporal sampling of the input SFHs. Our result is in 
general agreement with modeling of the spectral energy 
distributions of dwarf galaxies that show changes in the 
SFR on short timescales can account for d iscrepancies in 
the fl ux predictions of up to a factor of 2 (j.Tohnson et al.l 
1201,111 . The flux for one of the higher luminosity galaxies, 
NGC 6822, is consistently over-predicted by the models. 
This system is unique in the STARBIRDS sample as the 
star-formation activity has increased significantly over 
the most recent time bin of 20 Myr. It is possible that 
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Figure 9. Ratios of the predicted to observed NUV fluxes as a 
function of the B-band absolute magnitude. The predictions were 
generated using the CMD-based SFRs and the synthetic spectral 
model (top left), synthetic stellar populations (top right), a syn¬ 
thetic galaxy taking into account stochastic effects (bottom left), 
and a synthetic galaxy spectrum (bottom right). Starburst galax¬ 
ies are plotted as filled symbols; post-starburst galaxies are plotted 
with unfilled symbols. The dashed line represents the mean NUV 
ratio and the shaded grey regions denote the standard deviation on 
the mean. The galaxies hosting young, massive clusters are plotted 
as stars and are excluded from the calculation of the mean ratio. 

the SFR has fluctuated significantly during this shorter 
time period, which would affect the flux predictions. 

The overall agreement between the NUV flux predic¬ 
tions from the four different models and the measured 
fluxes provides an independent check on a number of vari¬ 
ables. First, the UV attenuation has been adequately 
corrected in the GALEX fluxes using the infrared flux 
from the MIPS images. If the measured fluxes were sig¬ 
nificantly affected by extinction at UV wavelengths, the 
predictions would be systematically biased high. Sec¬ 
ond, the SFHs based on GMD fitting are accurately mea¬ 
suring the SF at recent times in these systems. Third, 
the transformations from the stellar evolution libraries to 
NUV wavelengths are consistent with what is measured 
from the composite stellar populations. Fourth, stochas¬ 
tic sampling of the IMF at lower SFRs does not appear 
to be a significant factor in our results as the predic¬ 
tions from both Pegase and GALAXEV, which do not 
include stochastic effects, are in agreement with the mea¬ 
sured values and consistent with predictions from MATCH 
and SLUG, which do include include such effects. This is 
not surprising as the stellar populations responsible for 
producing the majority of the NUV fluxes includes the 
comparatively more abundant stars of lower masses than 
the FUV emission, reducing the impact of whether or not 
the upper end of the IMF was fully populated. 

4.2. FUV Ratios 

Figure fTOl shows a similar analysis to Figure El for the 
FUV bandpass. Here, the predicted FUV fluxes from 
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Ratios of Predicted to Measured FUV Fluxes 
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Figure 10. Ratios of the predicted to observed FUV fluxes as a 
function of the B-band absolute magnitude. The predictions were 
generated using the CMD-based SFRs and the synthetic spectral 
model (top left), synthetic stellar populations (top right), a syn¬ 
thetic galaxy taking into account stochastic effects (bottom left), 
and a synthetic galaxy spectrum (bottom right). Starburst galax¬ 
ies are plotted as filled symbols; post-starburst galaxies are plotted 
with unfllled symbols. The dashed line represents the mean FUV 
ratio and the shaded grey regions denote the standard deviation on 
the mean. The galaxies hosting young, massive clusters are plotted 
as stars and are excluded from the calculation of the mean ratio. 

the four models are seen to disagree with the measured 
GALEX FUV fluxes. The mean predicted to measured 
FUV ratio for Pegase is low with x = 0.54, cr = 0.22, 
whereas the mean FUV ratios for MATCH, SLUG, and 
GALAXEV are consistent with one another but high 
ranging from x = 1.54, cr = 0.17 to a; = 1.56, cr = 0.13 
and X = 1.74, CT = 0.28. The SLUG results show a lower 
dispersion than the other models. This is most likely 
due to the inclusion of additional stochastic effects when 
clustered star formation is taken into account. 

The offsets in the predicted FUV fluxes are likely due 
to the calibration of the models at the FUV bandpass, 
which depends on both the application of stellar atmo¬ 
sphere libraries and the calibration of the stellar evo¬ 
lutionary isochrones at FUV wavelengths. As noted in 
Section 13.2.51 there are important differences between 
the four models that impact the calibrations at FUV 
wavelengths. All four codes use different stellar atmo¬ 
sphere models in the UV regime as summarized in Ta- 
ble[5] Pegase, SLU ( j. GA LAXEV use the models from 
iLei^ne et al.l (|1997L 1199811 (based on the original grid s 
of model atmosph ere sp ectra from iBessell et al.l (1198911 . 
iFluks et"^ (1199411 , and iKuruczl (1199211 1 for the major¬ 
ity of stars contributing to the galaxy spectra. However, 
SLUG and GALAXEV use additional models for stars 
with strong stellar winds, O stars, and non-solar metallic- 
ity sta rs. MATCH uses the library from iGastelli &: Kurucz 
(2004), that has some overlap with the original iKuru^ 
( 199211 model spectra. In addition, MATCH uses a zero 
point transformation for the isochrones to the GALEX 


FUV bandpass. 

While disentangling the different effects of these li¬ 
braries and zero point calibrations is beyond the scope 
of this study, we did perform a series of simple tests to 
measure the impact stellar atmosphere models can have 
on FUV fluxes. We ran a series of simulations using 
the Starburst99 spectral models, which allows the user 
to model a galaxy spectrum with different atmospheric 
libraries. We used the same stellar evolutionary libraries 
and IMF as our FUV predictions, keeping the inputs and 
assumptions as closely matched as possible. The simu¬ 
lations spanned a SFR range from 0.001 to 1.0 Mq yr“^ 
and a metallicity range from Z of 0.0004 to 0.02. While 
the use of another code may introduce additional uncer¬ 
tainty in the comparison, the Starburst99 code provides 
a relative measure of how the stellar atmosphere models 
affect the FUV fluxes. Further, not only is the technique 
used in Starburst99 similar to GALAXEV, the SLUG 
code is based on the Starburst99 code, with modifica¬ 
tions to include stochastic effects. Thus, the additional 
uncertainties in the relative comparison should be small. 

Starburst99 simulations were run using the atmo¬ 
spheres applied in SLUG and the atmospheres used in 
Pegase with the m atched stellar evolution isochrones 
(|Bressan et al.l[T99^ : the atmospheres used in MATCH and 
GALAXEV were not available in the Sta rburst99 simu- 
latioii s . The Sta rburst99 fluxes using t he iLeieune et al.l 
(jl997l [T998fl and iHillier fc Miller! (jl998D models, similar 
to the SLUG stellar atmosphere models, varied by up to 
10 % at solar metallicities, with smaller differences at the 
lower metallicities of the STA RBIRDS samp le. The Star- 
burst99 fluxes using just the ILeieune et al'l (|1997L I1998D 
models with the older Padova models, similar to Pegase, 
varied by a similar amount with a similar dependency on 
metallicity. 

While our simple test using Starburst99 did not pro¬ 
duce big enough changes in the FUV fluxes to account 
for the discrepancies seen in Figure 1101 we cannot rule 
out differences in the application of the atmospheres as 
a sign ificant source of the offset. iBruzual fc Charlatl 
(120031) note in their Appendix that the UV colors of the 
spectra are significantly more uncertain than the optical 
colors due to the l ack of comparison st andards at non¬ 
solar metallicities (|Westera et al.l 1200^ . They qualita¬ 
tively conclude that the reliability of the predictions for 
young stars at low metallicity is fair to poor. A previous 
comparison between P egase and Starburs t99 reported 
significant differences (|Kewlev et al.l 1^00ll) . However, 
not only were the model atmospheres different, but the 
earlier version of the Starburst99 code us ed a different 
stellar evolutionary library. Thus, while iKewlev et al.l 
(j2001fl speculated that the different model atmospheres 
were driving the differences, the different evolutionary li¬ 
braries prohibited isolating the effects of just the model 
atmospheres. 

5. COMPARISON OF SFRS FROM INTEGRATED FUV 
EMISSION USING EXISTING SGALING RELATIONS 
WITH SFRS(T) FROM RESOLVED STELLAR 
POPULATIONS 

Figure [TT] presents a comparison of the SFRs measured 
from the resolved stellar populations averaged over the 
past 100 Myr to the SFRs based on the extinction cor¬ 
rected, integrated FUV emission using the scaling rela- 
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Figure 11. Comparison of the CMD-based SFRs averaged over 
the last 100 M yr with FUV-base d SFRs calculated using the scaling 
relation from IHao et al.l II2011J) . The uncertainties on the CMD- 
based SFRs are based on the standard deviation in SFRs averaged 
over 50, 100, and 150 Myr, thus providing a measure of how variable 
the SFRs have been over the timescale of interest. Starburst galax¬ 
ies are plotted as filled circles; post-starburst galaxies are plotted 
as unfilled circles. The solid black line indicates a ratio of unity 
between the two SFRs; the dotted line is where agree ment would 
be us ing the FUV-based SFRs scaling relation from ISalim et al.l 
H2007l'l . The CMD-based SFRs are higher than the FUV-based 
SFRs for all systems except in four cases, plotted as stars, whose 
CMD-based SFRs are lower-limits and which also compromise the 
assumption of constant SF at the most recent times as discussed 
in the text. 

tion from iHao et al.l (|2011f) . As noted above, the uncer¬ 
tainties on the CMD-based SFRs are the standard de¬ 
viations of the SFRs between timescales of 50, 100, and 
150 Myr. We included the time bin of 150 Myr in the 
standard deviation as it extends the temporal baseline 
over which assumption of constant SFR can be evaluated. 
The uncertainties in the FUV-based SFRs include statis¬ 
tical uncertainties only. For completeness, we include the 
four galaxies with young, massive clusters whose CMD- 
based SFRS are lower-limits (plotted as stars). As ex¬ 
pected, these galaxies have higher FUV-based SFRs than 
the CMD-based SFRs. 

The FUV-based SFRs are clearly correlated with 
CMD-based SFRs. However, the majority of the FUV- 
based SFRs are offset to lower values by up to a factor of 
4. The comparison to the FUV-based SFRs is affected by 
our choice of the scaling relation we used to convert from 
the FUV emis sion to a SFR Isee (13.41) . If we had used the 
relation from ISalim et al.l (|2007ll . the FUV-based SFRs 
would be 25% lower, represented by the dotted line in 
Figure 1111 This difference would increase the discrep¬ 
ancies between the CMD-based and FUV-SFRs. Fur¬ 
ther, in deriving the empirical s caling relation to co nvert 
the FUV luminosity to a SFR, ISalim et al.l (|2007f) cali¬ 
brated the results using SFRs derived from fitting spec¬ 
tral energy distributions from the UV to optical bands. 
In their work, the empirical SFR relation is noted to be 


a poorer fit to the SED-derived SFRs for values below 
^0.1 Mq yr“^; the empirical relation generally produc¬ 
ing FUV based SFRs that are higher than SED based UV 
SFRs with a significant amount of scatter. The authors 
attribute the differences to difficulties in accurately esti¬ 
mating the extinction based on FUV-NUV colors. Thus, 
not only do our results show lower FUV-based SFRs for 
the majority of the galaxies, but our analysis also cor¬ 
rects for extinction based on the far IR emission from 
dust in the galaxies. 

On the one hand, SFRs calculated from integrated 
FUV emission measurements are useful because of their 
simple application. On the other hand, interpreting the 
SFRs calculated in this way can be limited for the same 
reason. Without additional information (that can be dif¬ 
ficult to obtain for large surveys or distant galaxies), 
the sources of ambiguity in the interpretations include 
correcting for FUV attenuation from dust, evaluating 
whether the assumption of constant SFR is valid for a 
given system, and determining whether the data are fully 
sampling the IMF and/or the SFH of the system. Indeed, 
many of these variables are often invoked as plausible 
explanations for discrepancies in observational compari¬ 
son of integrated light SFRs. While these challenges are 
present in this study, the advantage of comparing inte¬ 
grated FUV flux measurements with SFHs derived from 
resolved stellar populations is that a number of the vari¬ 
ables can be more fully constrained. In Appendix A, we 
discuss several of these factors in detail and provide ar¬ 
guments to rule out these often invoked uncertainties as 
the source of the offset between the SFR indicators. 

The direction and amount of offset between the FUV- 
based and CMD-based SFRs is consistent with the offset 
between the FUV flux predictions from three previous 
models and the extinction corrected FUV emission mea¬ 
sured from the GALEX images (see Section 4.2). This 
is in contrast to the NUV flux predictions using the 
same inputs which were in very good agreement with 
the measured NUV emission (see Section 4.1). Two 
of the three models, GALAXEV and SLUG, have pre¬ 
viously been_^hov 2 ito_te consistent with Starburst99 
model (jda Silva et al.l[^12fl whic h is used to calibr ate the 
FUV-SFR scaling relation from iHao et al.l (1201111 . This 
strongly suggests that the most likely cause of the offset 
between both the two SFR indicators and the predicted 
to measured FUV fluxes is the calibration of the FUV 
flux used in the models, either from the rather uncertain 
stellar atmosphere models for UV producing stars, a sys¬ 
tematic in the stellar evolutionary isochrones at FUV 
wavelengths, or both. The lower FUV-based SFR for a 
given GMD-based SFR implies that the theoretical mod¬ 
els used in Starburst99 to calibrate the FUV-SFR scaling 
relation over-predict the FUV flux for a given SFR. 

6 . EMPIRICAL CALIBRATION OF THE FUV-BASED SFR 
RELATION USING CMD-BASED SFRS 

Here, we quantify the difference between the CMD- 
based SFRs and the FUV-based SFRs for the STAR- 
BIRDS sample. We present an alternative calibration 
for converting integrated FUV emission to a SFR us¬ 
ing the galaxies from the STARBIRDS sample that best 
meet the assumption of constant SF activity over the 
lifetime of the FUV emitting stellar population. First, 
we eliminate the four galaxies that host young stellar 
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Figure 12. The CMD-based SFRs are compared with the ex¬ 
tinction corrected FUV luminosities. The red circles represent the 
galaxies whose SF activity is fairly constant over 100 Myr; un¬ 
certainties in the HST SFRs are approximat ely t he size of the plot 
points. Individual errors are plotted in Figure flTI The least squares 
fit to these points is sho wn as a solid red lin e. For comparison, we 
plot the calibration fromISalim et al. 1120071) as a dash-dot line and 
the calibration from lHa'oe^dTl20n ) as a dotted line. The differ¬ 
ent values of the calibrations are noted in the Figure legend. The 
black stars represent galaxies with young, massive star clusters 
whose CMD-based SFRs are under-reported and which also com¬ 
promise the assumption of constant SF at the most recent times. 
The blue triangles represent the galaxies that show variations in 
the SF activity at older times. While these are excluded from the 
calibration, their SFH variations do not impact their FUV emis¬ 
sion as much as in systems with significant variations within the 
last 10 Myr. 

clusters with masses of order 10 ^ Mq that both com¬ 
promise the assumption of constant SF and have unre¬ 
solved stars in the clusters in the optical images (i.e., 
NGC 4449, NGC 5253, NGG 4214, and NGG 625). Sec¬ 
ond, we eliminate three galaxies that show significant 
variability in their GMD-based SFHs over the past 100 
Myr, namely Antlia, UGG 9128, and the newly bursting 
galaxy NGG 6822. 

We adopt th e notation for the FUV-base d SFR relation 
introduced by iGharlot fc Longhettil (120011) where the lu¬ 
minosity of a galaxy is related to a SFR by the efficiency 
factor rj: 

Vfuv ~ ^%uv/ (CMD-based SFR)iqq Myr (9) 

In Figure [121 we plot the CMD-based SFRs and the 
extinction corrected FUV luminosities for the sample. 
Using the galaxies that best approximate constant SFR 
over the last 100 Myr (red circles), a least squares fit 
including uncertainties on the CMD-based SFRs yields a 
value of log 'r]p^Y = —27.69. The inverse of rjpjjy gives 
the conversion factor for the FUV-based SFR as: 

SFi? = 2.04 ± 0.81 X 10-28 X (10) 

where SFR is in Mq yr”^ and Apc/v ^^e extinction cor¬ 


rected luminosity measured in the GALEX FUV band¬ 
pass in erg s”^ Hz”^. This calibration is based on a 
randomly populated Salpeter IMF with mass limits be¬ 
tween 0.1 and 120 Mq, for low-metallicity dwarfs with 
SFRs between 10“^ and 10“^ Mq yr”^. T his is 1.53 
times greater than the calibration reported in iHao et al.l 
(|201lD an d 1.72 times great er than the calibration re¬ 
ported in I Salim etldl (120071) for higher-metallicity and 
higher-SFR galaxies. Our empirically determined scal¬ 
ing factor includes an estimate of the uncertainty based 
on the standard deviation of the data from the best fit 
line, including uncertainties on the individual measure¬ 
ments. This offers an improvement over theoretical cali¬ 
brations that do not quantify the uncertainties. Table [3] 
lists the FUV SFRs for our sample calculated with the 
new calibration in Equation 1101 

Note that our sample is small, and therefore the cal¬ 
ibration presented here could be improved upon with a 
larger number of galaxies. However, even with a larger 
sample of galaxies, there will likely still be a higher level 
of dispersion in the calibration in the low-SFR regime 
of dwarf galaxies. To test whether the dispersion in the 
correlation between SFR to FUV luminosity is due to 
observational uncertainties or is intrinsic to the relation¬ 
ship, we fit a linear relationsh ip using the IDL mpfitexy 
routine (|Williams et al.ll201Cll see Section 4). The mpfi¬ 
texy routine depends on the mpfit package iMarkward^ 
(j2009f ). This procedure hnds a high degree of intrinsic 
scatter (97%) or dinstrinic = 0.79 x IQ-^® from Equa¬ 
tion [TOl The high intrinsic dispersion reflects the com¬ 
plex nature of these systems, but also note that this value 
would be lower if the reported uncertainties are under¬ 
estimated. We discuss the dispersion between the SFR 
and FUV luminosity further in the next section. 

Similar finding s _ have been reported by 

iGhomiuk &: Povichl (|2011l) in a comparison of SFRs mea¬ 
sured in the Milky Way with a variety of SFR indicators. 
These authors report that Galactic SFRs derived for 
H II regions using resolved low and intermediate mass 
stellar populations were systematically higher by factors 
of ^ 2 to 3 than SFRs derived using high mass stellar 
tracers in the mid-IR and radio emission. A separate 
study reported SFRs derived from H/3 line strengths to 
be ^ 1.8 X higher than those derived fro m UV fluxes in 
a sam ple of galaxies at a redshift of ~ 2 (jZeimann et al.l 

(Ml) . 

In contrast, comparisons of integrated FUV-based 
SFRs to SFRs derived from SED fits ha ve found the 
FUV- based SFRs to yield higher SFRs. iHuang et al.l 
(j2012D compared integrated FUV-based SFRs with SED 
fits of a sub-sample of low-mass galaxies from the Arecibo 
Legacy Fast ALFA (ALFALFA) Survey. These authors 
report that the FUV-based SFRs were slightly higher 
than the SFRs from SED htting, in general agreement 
with the results fr om a larger s ampl e of more massive 
galaxies studied bv I Salim et al.l (j2n07[ ). 

7. THE NATURE OE THE DISPERSION IN FUV-BASED 

SFRS 

The CMD-based SFRs and extinction corrected FUV 
luminosities shown in Figure [T^] are clearly correlated 
but there is still a high degree of dispersion. The un¬ 
certainties in the CMD-based SFRs and the statistical 
uncertainties in the FUV luminosities, which are approx- 
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Figure 13. Top left: Histogram of the ratio of CMD-based to 
FUV-based SFRs using the new calibration. Top right: The CMD- 
based SFRs versus the ratio of the CMD-based SFRs to FUV-based 
SFRs. The points are color-coded by metallicity as described in the 
legend. The is no correlation seen between the dispersion in the 
SFR indicators and the SFR. Bottom left: The percent of diffuse 
FUV emission versus the ratio of CMD-based SFRs to FUV-based 
SFRs. Uncertainties represent the difference in diffuse emission 
if the empirically selected threshold delineating diffuse from UV- 
bright knots is varied by ±50%. There are no galaxies with a high 
fraction of diffuse emission with lower CMD-based to FUV-SFR 
ratios or vice-versa. Bottom right: Maximum surface brightness of 
the galaxies versus the ratio of CMD-based SFRs to FUV-based 
SFRs. Similar to the previous panel, there are no galaxies with 
higher surface brightness values and high CMD-based to FUV- 
based SFRs, or vice-versa. 

imately the size of the points, generally do not overlap 
with the best-fit line used for the calibration. We plot 
the differences between the CMD-based SFRs with the 
FUV-based SFRs calculated using our new calibration 
in the upper left panel of Figure ITSl For clarity, we plot 
only the galaxies used in the SFR calibration (see Sec¬ 
tion 6). The ratios of the two SFR indicators reveal that 
there is still up to a factor of two difference between the 
measurements. 

As discussed in Section 6, most of this dispersion ap¬ 
pears to be intrinsic and can be due to a number of differ¬ 
ent factors. The more likely factors include (i) stochastic 
effects on the upper-end of the IMF which depends on 
the level of star-formation activity, (ii) the known depen¬ 
dency of the FUV flux on stellar metallicity, and (iii) the 
amount of concentrated versus diffuse FUV emission. In 
Figure [131 investigate the dependency of the disper¬ 
sion in the SFRs on these three factors using measurable 
quantities, namely the level of star-formation activity, 
gas-phase metallicity measurements, and the amount of 
concentrated versus diffuse FUV emission. 

In the top right panel of Figure [T31 we plot the ratio 
of CMD-based SFR to FUV-SFR calculated using our 
new calibration constant against the CMD-based SFRs. 
Based on the premise that higher SFRs will be less sus¬ 
ceptible to stochastic effects, the expectation is that there 
will be an decrease in dispersion at higher SFRs. How¬ 
ever, there is no clear correlation between the dispersion 
and the star-formation activity, suggesting that this is 
not the dominate driver of the dispersion. Note that this 


does not imply that stochastic effects do not contribute 
to the scatter. As demonstrated in Section 4.2, the FUV 
flux predictions from the MATCH and SLUG models — 
which stochastically populate the IMF — have the low¬ 
est scatter of the models tested. 

Second, we consider the known dependency the metal¬ 
licity of a star has on the FUV emission. At a given 
SFR, there is less UV flux produced from higher metal¬ 
licity stellar populations due to lower stellar tempera¬ 
tures and line blanketing. Thus, when a constant scal¬ 
ing relation is applied over a range in metallicity, there 
is an increasing deviation between the input and recov¬ 
ered SFRs as a function of metallicity. Previous authors 
have reported on the amplitude of this dependency, with 
UV fluxes varying by less than 0.24 dex, or ^ 70%, 
ove r a metallicity range from Z of 0.0003 to 0.03 (e.g., 
see IMadau fc PickinsM 1201411 . Other authors report a 
change o f 55% in the flux o ver log Zq /Z range of -1-0.2 
to —1.0 (iConrov et al.ll200'^ . 

Despite this known dependency, we do not find an em¬ 
pirical signature for the FUV-based SFR calibration to 
be changing as a function of metallicity. In Figure [T51 
the points are color-coded by metallicity bins with bluer 
colors representing galaxies that are more metal-poor 
and yellow colors representing galaxies that are compara¬ 
tively more metal-rich (see legend for metallicity ranges). 
Changes in the UV flux as a function of metal content 
over this range are likely present in the stellar popula¬ 
tions, yet the dispersion between the SFR indicators do 
not show a clear dependence on metallicity. From our 
Starburst99 simulations (see Section 4.2), the expected 
amplitude of the dependence is ^ 15% over the narrow 
metallicity range of the STARBIRDS sample. This de¬ 
pendency is smaller than the offset between the CMD- 
based SFRs and the FUV-based SFRs and, thus, is prob¬ 
ably not the dominant driver of the dispersion. 

Third, inspection of the GALEX and HST images 
of the galaxies suggests that the presence and num¬ 
ber of star forming regions with concentrated UV emis¬ 
sion (“UV-bright knots”) and the fraction of emission 
originating outside of such regions (“diffuse emission”) 
may be more directly correlated with the dispersion 
in the FUV-based SFRs. The UV-bright knots trace 
the young star-forming complexes in the galaxies; the 
source of diffuse UV emission is not clear. On the one 
hand, some of the diffuse UV light may originate from 
the main star forming complexes and be scattered by 
dust and gas within a galaxy. On the other hand, re¬ 
cent studies have suggested that the diffuse emission is 
from B stars that have migrated from stell ar associa¬ 
tions into field regions of the galaxies (e.g.. iCole et al.l 
119991 : iTremonti et al.ll2001t ICliandar et al.l I^OSfl The 
timescale required by these studies for dispersing the 
B stars from the star forming associations is of order 
^ 10 Myr. However, other studies report a substan¬ 
tially longer timescale of ^ 200 Myr for st ellar clusters 
in dwarf galaxies to dissolve in to the field (iGieles et al.l 
120081 iBastian et ^I20ini2011h . This result is supported 
by a separa te study on the cluste r disruption f raction in 
NGC 4449 (jAnnibali et al.ll2011h . Separately. iLiu et al.l 
(|2011ll report that 80% of the FUV emission originates 
from stars with ages < 30 Myr compared to 20% of the 
FUV emission from older, lower-mass stars that have 
presumably migrated from their clustered birth places. 
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However, the actual measured fraction of emission in 
clusters can be as low as 60% due to the higher levels 
of attenuation in the more dusty nascent regions . This 
can be compared with work from I.Tohnson et al.l 1)20131 1 
who suggest that the mean age of the stellar population 
contributing to the FUV emission is strongly dependent 
on the SFH of a system, ranging from ~ 50% produced 
by stars younger than 16 Myr to 100% of the flux pro¬ 
duced by stars older than 100 Myr. 

Regardless of the origin of the diffuse emission, we in¬ 
vestigated whether the dispersion in the ratio of SFR in¬ 
dicators is due to the relative contribution of UV-bright 
knots to diffuse emission. We quantify the overall dif¬ 
ferences in the UV emission in two different ways from 
the GALEX data. First, we consider how the fraction 
of diffuse UV emission relative to the UV-bright knots in 
the STARBIRDS sample correlates with the difference in 
FUV-based SFRs. This is an approach that is only pos¬ 
sible for spatially resolved galaxies. Second, we compare 
how the maximum intensity in FUV emission correlates 
with the same difference, which enables an evaluation of 
this property for unresolved galaxies. 

In the bottom left panel of Figure [131 we plot the per¬ 
centage of diffuse emission compared to the ratio of the 
CMD-based SFRs to FUV-based SFRs. A rough esti¬ 
mate of the fraction of diffuse emission was calculated 
from the ratio of the flux in all pixels less than an em¬ 
pirically selected threshold of 0.025 counts s“^ pixeU^ 
to the total flux in the FUV image. The threshold value 
was determined by examining the flux contour levels that 
encompass the UV-bright knots across the sample. The 
uncertainties represent the range in the percent of dif¬ 
fuse emission when the threshold value of 0.025 counts 
s“^ pixeU^ is varied by ±50%. As the UV-bright knots 
are sharply peaked in flux, varying this threshold does 
not significantly change the distribution of points. 

There appears to be a loose correlation in Figure |T2| 
which shows that galaxies with larger amounts of diffuse 
emission have GMD-based SFRs that are larger than the 
FUV-based SFRs. On one hand, this loose correlation 
could be the result of changes in the SFR at slightly 
older times compared to the most recent 100 Myr. On 
the other hand, the origin of diffuse FUV light is not 
well understood. A detailed comparison study of the re¬ 
solved stellar populations in both the UV and the optical 
might help clarify the source of the emission. However, a 
Kolmogorov-Smirnov test rejects the null hypothesis that 
these two samples are drawn from the same distribution 
with a p-value of < 1 %, indicating that these variables 
are not likely dependent upon one another. 

Because it is only possible to calculate the fraction of 
diffuse emission in resolved galaxies, we probed whether a 
different tracer of UV bright knots, the maximum inten¬ 
sity in FUV emission, could be correlated with the dis¬ 
persion in FUV-based SFRs. In the bottom right panel 
of Figure [T51 we plot the maximum value of counts s“^ 
pc“^ for each of the GALEX FUV images versus the ra¬ 
tio of GMD-based SFRs with FUV-based SFRs. While 
there is no clear trend in the data, there are no galaxies 
with higher values of both maximum surface brightness 
and GMD-based to FUV-SFR ratios or vice-versa. How¬ 
ever, similar to the previous comparison, a Kolmogorov- 
Smirnov test returns a a p-value of < 1% indicating these 
samples are likely drawn from different distributions. 


8 . CONCLUSIONS 

We have compared the star formation traced by extinc¬ 
tion corrected GALEX UV imaging in 19 resolved star- 
burst and post-starburst dwarf galaxies with HST optical 
imaging of the resolved stellar populations using data ac¬ 
cessible in the Panchromatic ST ARBIRDS Archive an d 
described in a companion paper (jMcQuinn et al.ll20I5fl . 
The UV emission was corrected for extinction by esti¬ 
mating the dust content from MIPS imaging. Here we 
summarize our conclusions. 

1. We find that SFHs derived from resolved stellar 
populations can accurately predict the integrated 
NUV fluxes in the GALEX bandpass using four in¬ 
dependent models for the majority of the sample. 
This result provides independent confirmation of 
the transformation of recent SFRs to NUV fluxes 
from stellar evolution libraries, the robustness the 
SFHs based on GMD fitting methods, and the ap¬ 
plication of UV extinction corrections in our sam¬ 
ple. In contrast, the integrated FUV fluxes pre¬ 
dicted using the same inputs and models do not 
agree with the extinction correction FUV fluxes 
measured from the GALEX images. The offset 
at FUV wavelengths is likely driven by calibra¬ 
tion differences in the stellar atmospheric models, 
systematic uncertainties in the stellar evolutionary 
isochrones, or both. 

2. For the majority of the sample, the FUV-based 
SFRs based on existing scaling relations are sys¬ 
tematically lower than the GMD-based SFRs. We 
show that the offset in the FUV-based SFRs can¬ 
not be explained by uncertainties of variable SFHs, 
mismatched emission timescales, or attenuation by 
dust. Additionally, while stochastic sampling of 
the IMF at low SFRs, particularly in highly clus¬ 
tered SF, introduces dispersion in the SFRs cal¬ 
culated using FUV-based SFR scaling relations in 
this SFR regime, stochastic effects cannot account 
for the discrepancy between the FUV-based SFRs 
from GMD-based SFRs. We caution that our sam¬ 
ple is small; our results would benefit from further 
work on a larger sample of galaxies. 

The offset between the GMD-based SFRs and inte¬ 
grated FUV-based SFRs is in agreement with the 
offset between predicted to measured FUV fluxes 
using three different models. The offset is likely 
driven by the theoretical calibrations of either the 
stellar atmosphere libraries or stellar evolutionary 
isochrones (or both) used in Starburst99 to cali¬ 
brate the FUV-based SFR scaling relation and in 
the codes GALAXEV, SLUG, and MATCH used to 
predict the FUV fluxes. 

3. Using the GMD-based SFRs from the low- 
metallicity galaxies that best meet the assumption 
of constant star formation over the past ^ 100 Myr, 
we calculate a conversion from integrated GALEX 
EUV emission to EUV-based SFRs of 2.04 ±0.81 x 
10-28 This is the first empirical calibration of 
the FUV-based SFR scaling relation using a sam¬ 
ple with well-measured and well-constrained star- 
formation properties outside of the synthesis model 
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approach. The value of the conversio n factor is 
53% larger than that of lHao et al.l (1201 If) calibrated 
from the Starburst99 model. 

Because CMD-based SFRs are derived using more 
information from the stellar populations and thus 
use fewer assumptions, they are likely a truer mea¬ 
sure of the star-formation activity than the FUV- 
based SFRs calculated using theoretically cali¬ 
brated scaling relations. Additionally, the CMD- 
based SFRs use a stellar evolutionary library with 
transformations to the HST filter system. These 
transformations to optical wavelength use spec- 
trophotometric standards and lar ge libraries of 
stellar spectra (iGirardi et al.l 120081 and references 
therein). In contrast, the theoretical model used 
to calibrate the FUV-based SFR scaling relation 
relies on model atmospheres from a smaller num¬ 
ber of massive stars that are not as well understood 
as the lower-mass stars included in the HST trans¬ 
formations. Further, in contrast to the theoreti¬ 
cal calibrations, the empirical derivation allows a 
robust estimate of the uncertainty on the scaling 
factor. Thus, our empirically derived FUV-based 
SFR scaling relation promises an improvement over 
previous theoretical relations. 

4. The majority of the scatter between CMD-based 
SFRs and FUV emission from low-mass galaxies is 
intrinsic, reflecting the complex nature of these sys¬ 
tems. While stochastic sampling of the upper-end 
of the IMF likely introduces some scatter into the 
relationship, we do not find empirical signatures 
of a dependency on star-formation activity driving 
the dispersion. Further, we do not find an empirical 
signature of the known dependence of FUV emis¬ 


sion on metallicity driving the dispersion. We do 
note a possible dependency of the conversion factor 
on the fraction of diffuse FUV emission compared 
to UV-bright knots in a galaxy, although statistical 
tests showed little probability of a correlation. A 
detailed understanding of the source of the diffuse 
emission is needed before its impact on measuring 
SFRs from integrated light can be further quanti¬ 
fied. 

9. ACKNOWLEDGMENTS 

Support for this work was provided by NASA through 
an ADAP grant (NNX10AD57G), and a NASA GALEX 
grant (NNX10AH98G). N. P. M is grateful for an NSF 
REU grant (PHYS 0851820) at the University of Min¬ 
nesota. GALEX (Galaxy Evolution Explorer) is a 
NASA Small Explorer, launched in 2003 April. We 
gratefully acknowledge NASA’s support for construction, 
operation, and science analysis for the GALEX mis¬ 
sion, developed in cooperation with the Centre National 
d’Etudes Spatiales of France and the Korean Ministry 
of Science and Technology. This research made use of 
NASA’s Astrophysical Data System, the NASA/IPAC 
Extragalactic Database (NED) which is operated by the 
Jet Propulsion Laboratory, California Institute of Tech¬ 
nology, under contract with the National Aeronautics 
and Space Administration, and the HyperLeda database 
(http://leda.univ-lyonl.fr). Finally, we would like to 
thank Michel Fioc, Michele Fumagalli, Cai-Na Hao, 
Claus Leitherer, and Brigitte Rocca-Vohnerange for help¬ 
ful discussions on modeling star formation, and John 
Cannon and Yaron Teich for constructive feedback on 
the manuscript. 

Facilities: GALEX Space Telescope, Hubble Space 
Telescope, Spitzer Space Telescope 


APPENDIX 

POSSIBLE BUT UNLIKELY SOURCES OF THE OFFSETS BETWEEN FUV-BASED SFRS AND CMD-BASED SFRS 

Comparisons between SFR indicators can be difficult because each indicator is based on a variety of assumptions 
and can be impacted by different observational factors. For example, an FVU SFR is based on the predication that the 
SFR in a system has been constant over the lifetime of the FUV emitting stars. CMD-based SFRs assume an accurate 
accounting of stars above the completeness with minimal loss in the photometry due to blending and source confusion. 
Differences in attenuation in the wavelength dependent measurements introduce uncertainty in SFR comparisons and 
are often invoked as an explanation when optical and UV SFR indicators do not fully agree. These make up a few of 
the possible sources of the offset between the FUV-based SFRs and the CMD-based SFRs reported here. However, the 
advantage of the STARBIRDS data and this study is that many of these assumptions and differences can be explored 
in detail. Here, we itemize a number of possible explanations for the discrepancy between the FUV and CMD-based 
SFRs and discuss why these are unlikely sources of the offset seen in Figure ITT] 

1. The deviations in FUV-based SFRs cannot be attributed to extinction correction uncertainties at the UV wave¬ 
lengths. This is supported by several arguments. First, UV attenuation is expected to be low as the sample 
is mainly comprised of low metallicity dwarfs. Second, uncertainties due to differences in metallicity and UV 
extinction would more likely introduce scatter into the relationship, instead of an offset. Third, there is no corre¬ 
lation between the ratio of the CMD-based SFRs to FUV-based SFRs and the oxygen abundance of the galaxies. 
Fourth, careful accounting of the UV extinction was performed using measurements of the infrared emission from 
the MIPS images. Finally and most importantly, the FUV fluxes predicted from the resolved stellar populations 
were in good agreement with the measured, extinction corrected FUV fluxes. Since the CMD-based SFRs are 
not as heavily impacted by extinction as the UV fluxes, the agreement between predicted and measured extinc¬ 
tion corrected UV fluxes independently confirms that the UV fluxes used in the SFR scaling relations are not 
significantly impacted by extinction. 

For the four galaxies whose FUV-based SFRs are higher than the CMD-based SFRs, we consider the possibility 
that we have oi;er-corrected for extinction. This could occur if the assumption that the majority of the UV 
flux comes from young stellar populations does not hold. If the primary source for heating the dust originates 
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from other mechanisms, such as mechanical energy from stellar winds, galactic winds, and supernovae and/or 
light from somewhat older stellar populations, then using the IR emission as a measure of UV extinction will 
overestimate the correction. However, the extinction corrections for these low metallicity dwarfs are fairly modest. 
Only by removing all of the extinction corrections from the GALEX fluxes do the measured values agree with 
the predictions. In other words, for the FUV-based SFRs to agree with the CMD-based SFRs, all of the IR 
emission would have to be attributed to mechanisms other than re-processed UV light. As this would be an 
unphysical interpretation of the IR emission, we rule out over-correcting the GALEX fluxes for extinction as the 
main source of the discrepancy in these four systems. 

2. The deviations in the FUV-based SFRs from the CMD-based SFRs are not dominated by under-sampling the 
upper-end of the IMF at low SFRs. As discussed above, while SFR scaling relations are based on the assumption 
of a fully-populated IMF, the upper-end of an IMF may not be fully populat ed at low-SFRs. In the UV continuum, 
this introduces an uncertainty of only ^20% for SFRs ^ 0.0003 M 0 yr’ (|Lee et alJl20Tni . Further, if the lower 
FUV-based SFRs were a result of less high-mass stars at lower SFRs depressing the expected UV flux, there 
would a dependency in the offset on the SFR. The SFRs shown in Figure [TT] span nearly 4 orders of magnitude, 
but the offset is roughly constant across the whole sample (excluding the galaxies with massive clusters previously 
discussed in Section 3.1.1). Thus, while stochast i city does have an effect on the measured UV fluxes as de scribed 
by numerous studies (e.g., iGoddard et al.l[2f)Tol : iLee et all 120111: iKoda et'n]l 2012 t l.lohnson et al.l[20T^ and is 
likely a contributing factor to the scatter in the comparison, it does not explain the offset between SFRs measured 
in the sample. This is supported by our comparison between the predicted fluxes based on the CMD-based SFRs 
and the three separate models discussed in Section 4. In contrast to the MATCH and SLUG models which randomly 
populate an IMF for a given SFR, the GALAXEV model does not. Yet all three are in agreement with one 
another. 

3. The offset in the FUV-based SFRs is not due to a breakdown in the assumption of constant star formation. 
Scaling the FUV emission to a SFR is sensitive to the assumption of constant star formation activity over the 
~ 100 Myr lifetime of the FUV emitting population of stars. When considering the uncertainties, this assumption 
of constant star formation is generally supported in the STARBIRDS sa mple. We can eva l uate th is assumption 
on a galaxy by galaxy basis. Looking at the SFHs shown in Figure 2 of lMcQuinn et aO (|2010bl l. a number of 
galaxies have only small variations in their SFRs over the past 100 Myr. Specifically, NGC 4068, ESO 154—023, 
NGC 784, DDO 165, NGC 6789, NGC 4163, and NGC 1569 have experienced sustained star formation activity 
over this timescale, yet show the largest discrepancy with the FUV-based SFRs of more than a factor of 4 below 
the CMD-based SFRs. A few of the galaxies, such as UGC 4483, UGC 6456, NGC 2366, and Holmberg II, do 
deviate from the assumption of constant star formation over the past 100 Myr, with increases in their SFR within 
the most recent ^ 20 — 40 Myr. The FUV-based SFRs for these galaxies are not as low relative to the CMD-based 
SFRs, reflecting this recent increase in SF activity. The deviation from a constant SFR actually improves the 
agreement between the SFRs. Note that while the CMD-based SFHs track changes in star formation over the 
past 100 Myr, as with any technique, the temporal sampling achievable is finite. One can imagine changes in the 
SFRs on shorter timescales than are recovered in the SFHs. As discussed above, such changes do occur and will 
introduce scatter into the comparison between the integrated FUV emission SFRs and the CMD-based SFRs. 
However, contrary to these fluctuations that are expected to introduce scatter to the relationship, we find an 
offset. 

4. The offset between the FUV-based SFRs from the CMD-based SFRs are not due to differences in timescales of 
emission between the two wavelength regimes. There is still some debate on the timescale over which the FUV 
emission originates. For more ma ssive galaxies, models suggest tha t 90% of the FUV emission originates from 
stars less than 100 Myr old le.g.. iKennicuttl [19981 iHao et al.l 1201 111 . However, for low-mass galaxies there can 
be a significant contr ibution from stars older than 100 Myr to the integrated FUV flux, depending on the SFH 
(j.lohnson et al.l[20T3[) . Thus, we tested the sensitivity of the offset in FUV-based SFRs to different timescales. 
Specifically, we averaged the CMD-based SFR over timescales ranging from 50 - 200 Myr in increments of 50 
Myr. Within the timescales probed, the offset in FUV-based SFRs was present in the same systems with only 
slight variations. This can be seen in Figure [TT] where the uncertainties on the CMD-based SFRs are based on 
the standard deviation in the SFR fr om 50 — 1 5 0 My r. Very few of the error bars overlap with the FUV-based 
SFRs using the scaling relation from IHao et al.l (1201 It) . 

5. There was no dependency in the offset in the FUV-based SFRs with inclination angle. The geometry and 
inclination of a galaxy may impact the comparison for individual systems. For example, FUV-based SFRs 
derived from scaling relations in edge-on galaxies, such as NGC 784 and ESO 145—023, are known to have higher 
uncertainties. However, the UV fluxes predicted from the CMD-based SFRs were in good agreement with the 
measured, extinction corrected fluxes. Thus, while it may be that both the CMD-based and FUV-based SFRs 
are affected by the geometry in these galaxies, the inclination angle cannot alone explain the deviation between 
the SFRs. 
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Table 1 

Galaxy Sample and Properties 


Galaxy 

(1) 

R.A. 

(J2000) 

(2) 

Decl. 

(J2000) 

(3) 

Ms 

(mag) 

(4) 

Distance 
(Mpc) (i 

(5) 

D 25 

arcmin) 

(6) 

Ah 

(mag) 

(7) 

Oxygen 

Abundance 

(8) 

Ref. 

(9) 



Starburst 

Galaxies 





UGC 4483 

08:37:03.0s 

-|-69:46:31s 

-12.73 

3.2 

1.2 

0.091 

7.50 

1 

UGC 6456 

11:27:59.9s 

-|-78:59:39s 

-14.03 

4.3 

1.5 

0.096 

7.64 

2 

DDO 165 

13:06:24.85s 

-|-67:42:25s 

-15.09 

4.6 

3.5 

0.065 

7.80 

3 

IC 4662 

17:47:08.8s 

-64:38:30s 

-15.13 

2.4 

2.8 

0.188 

8.17 

4 

NGC 6822 

19:44:56.6s 

-14:47:21s 

-15.22 

0.5 

15.5 

0.632 

8.11 

5 

NGC 4068 

12:04:00.8s 

-|-52:35:18s 

-15.48 

4.3 

3.2 

0.058 

7.84 

3 

NGC 2366 

07:28:54.6s 

-|-69:12:57s 

-16.02 

3.2 

7.3 

0.097 

8.19 

6 

ESO 154-023 

02:56:50.38s 

-54:34:17s 

-16.40 

5.8 

8.2 

0.045 

8.01 

3 

NGC 784 

02:01:17.0s 

-|-28:50:15s 

-16.59 

5.2 

6.6 

0.158 

8.05 

3 

Ho II 

08:19:04.98s 

-|-70:43:12s 

-16.72 

3.4 

7.9 

0.086 

7.92 

7 

NGC 4214 

12:15:39.2s 

-|-36:19:37s 

-17.19 

2.7 

8.5 

0.058 

8.38 

8 

NGC 5253 

13:39:55.9s 

-31:38:24s 

-17.38 

3.5 

5.0 

0.186 

8.10 

9 

NGC 1569 

04:30:49.0s 

-|-64:50:53s 

-18.17 

3.4 

3.6 

1.871 

8.19 

10 

NGC 4449 

12:28:11.9s 

-|-44:05:40s 

-18.27 

4.2 

6.2 

0.051 

8.32 

11 



Post- 

■Starburst Galaxies 





Antlia Dwarf 

10:04:04.1s 

-27:19:52s 

-9.75 

1.3 

2.0 

0.212 

7.39 

12 

UGC 9128 

14:15:56.5s 

-|-23:03:19s 

-12.71 

2.2 

1.7 

0.065 

7.74 

13 

NGC 4163 

12:12:09.1s 

-|-36:10:09s 

-14.24 

3.0 

1.9 

0.052 

7.69 

3 

NGC 6789 

19:16:41.1s 

-|-63:58:24s 

-14.32 

3.6 

1.4 

0.187 

7.77 

3 

NGC 625 

01:35:04.6s 

-41:26:10s 

-16.53 

3.9 

6.4 

0.044 

8.10 

14 


References, — (1) [Skilljimn_et_a^ 


_ _ (5) |lriidalgo^^^tme^^^^m2uCTn ) 

( 10 ) [ 


j(ito3), (i4)isin! 


(2) Uroxal^i_ 

('6'> Ro^ (l'\ 

fc Ailfma Jliferr (11) 


ckv ik. Sk 

UIMM); 


Lee et al.l I2UU3 

): (8) iKobulnickv SkillmanI liuyfc 


Skillman et al.l 

19891); (12) Ipiersimoni et al.l 119991 

(13) 


Note. — Column (1) Galaxies are listed according to luminosity, from faintest to brightest. Columns (2) and (3) J2000 

coordinates. Column (4) luminosity corrected forextinctioi^JKara^ientse^^_e^ all2004ll .Column (5) Distance in Mpc. Column 

(6) Major axis of 25 mag isophote in arcmin (KaracEents^^^^ aljlsooj) . Column (7) Ajj (R=650 nm) extinction estimates 

are from the HI maps of [Schlege^_et_alJ (^993). Column (8) and (9) Oxygen abundance and reference. 
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Table 2 

Measured UV and IR Fluxes 


Galaxy 

(1) 

Fluxpuv X 10 Fluxjsjuv X 10 

(erg cm“^A“^) (erg cm“^A“^) 

(2) (3) 

Flux24/^m 

(Jy) 

(4) 

Flux7o/4m 

(Jy) 

(5) 

Fluxi6o/im 

(Jy) 

(6) 



Starburst 

Galaxies 




UGC 4483 


1.7±0.5 

0.87±0.13 

0.07 

0.1 

ND 

UGC 6456 


3.0±0.6 

1 .6±0.2 

0.03 


ND 

DDO 165 


6.4±0.9 

3.8±0.3 

ND 

0.06 

ND 

IG 4662 



16±1 

>5.0 

> 9.0 


NGC 6822:FOV 1 


3.7±0.7 

1.7±0.2 

0.2 

3.2 

6.0 

NGC 6822:FOV 2 


4.4±0.8 

2.3±0.2 

0.01 

3.2 

6.9 

NGC 6822:FOV 3 


5.6±0.9 

2.3±0.2 

0.2 

3.2 

4.0 

NGC 4068 


7.5±1.0 

4.1±0.3 

0.02 

0.7 

1.0 

NGC 2366 


31±2 

17± 1 

0.7 

5.0 

3.7 

ESO 154-023 


9.1±1.1 

5.1±0.3 

0.02 

0.7 

1.1 

NGC 784 


6.8±0.9 

4.7± 0.3 

0.03 

1.1 

1.4 

Ho II 


32±2 

17±1 

0.2 

2.7 

3.0 

NGC 4214 


54±3 

33±1 

1.7 

17 

19 

NGC 5253 


33±2 

24±1 

8.8 

26 

20 . 

NGC 1569 


3.2±0.7 

2 .8±0.2 

7.4 



NGC 4449 


160±5 

85±1 

3.1 

41 

71 



Post-Starburst Galaxies 




Antlia Dwarf 


0.33±0.22 

0.39±0.09 

ND 

ND 

ND 

UGC 9128 


1.2±0.4 

0.82±0.13 

ND 

ND 

ND 

NGC 4163 


2.5±0.5 

1 .8±0.2 

0.003 

0.07 

0.2 

NGC 6789 


0.8±0.3 

0.64±0.11 

0.008 



NGC 625 


9.7±1.2 

6.2±0.4 

0.8 

5.6 

5.8 


Note. — Column (1) Galaxy name. Columns (2)-(3) Integrated FUV and NUV flux measurements from the 
GALEX images in the HST fleld of view with foreground and background contamination masked. Uncertainties 
assume a Poisson distribution and are based on the square root of the number of counts at each wave band. 
Columns (4)-(6) MIPS flux measurements in the HST fleld of view after any evident foreground and background 
contamination has been masked. ND is a non-detection at that wavelength. Lower limits on the flux are noted 
for IC 4662 where the MIPS observations did not cover the HST field of view. 






In? 

to 


Table 3 

Measured UV Luminosities, Extinction Corrections, and Calculated FUV SFRs 


Galaxy 

(1) 

log Lpuv 
(ergs s“^) 
(2) 

log Lnuv 
( ergs s“^) 

(3) 

FUV 

(mag) 

(4) 

NUV 

(mag) 

(5) 

FUV-NUV 

(mag) 

(6) 

Apuv 

(mag) 

(7) 

Anuv 

(mag) 

(8) 

log 

(ergs s“^) 

(9) 

log ^%uv 

(ergs s“^) 
(10) 

log FUV SFR 
(M© yr-i) 
(11) 

log CMD SFR 
(M© yr-i) 
(12) 

log FUV SFR 
(M© yr-i) 
(13) 





Starburst Galaxies 








UGC 4483 

40.51 

40.38 

16.09±0.31 

16.02±0.17 

0.07 

0 

0 

40.51 

40.38 

-2.66 

-2.43 ± 0.02 

-2.47 ±0.21 

UGC 6456 

41.01 

40.90 

15.49±0.23 

15.36±0.12 

0.13 

>0 

>0 

41.01 

40.90 

-2.16 

-2.07 ±0.02 

-1.98 ±0.20 

DDO 165 

41.39 

41.33 

14.67±0.16 

14.43±0.08 

0.23 

0 

0.25 

41.39 

41.43 

-1.77 

-1.20 ±0.02 

-1.59 ± 0.18 

IC 4662 


41.39 


12.87±0.04 


>0.50 

>0.85 


41.73 


-1.52 ±0.01 


NGC 6822 (1) 

39.22 

39.05 

15.28±0.21 

15.31±0.12 

-0.03 

1.54 

1.38 

39.84 

39.61 

-3.33 

-2.86 ±0.04 

-3.14 ±0.17 

NGC 6822 (2) 

39.30 

39.18 

15.07±0.19 

14.98±0.10 

0.09 

1.41 

1.17 

39.87 

39.65 

-3.30 

-2.80 ±0.05 

-3.11 ±0.17 

NGC 6822 (3) 

39.41 

39.20 

14.82±0.17 

14.96±0.10 

-0.14 

1.11 

1.02 

39.85 

39.60 

-3.32 

-3.36 ±0.05 

-3.13 ±0.17 

NGC 4068 

41.40 

41.32 

14.50±0.15 

14.32±0.08 

0.18 

0.12 

0.02 

41.62 

41.33 

-1.72 

-1.47 ±0.02 

-1.53 ±0.18 

NGC 2366 

41.77 

41.67 

12.95±0.07 

12.80±0.04 

0.15 

0.40 

0.24 

41.92 

41.76 

-1.25 

-1.30 ±0.02 

-1.06 ±0.17 

ESO 154-023 

41.75 

41.66 

14.29±0.13 

14.12±0.07 

0.18 

0.07 

0 

41.77 

41.66 

-1.40 

-1.11 ± 0.01 

-1.21 ±0.18 

NGC 784 

41.53 

41.54 

14.61±0.16 

14.18±0.07 

0.43 

0.36 

0.11 

41.67 

41.59 

-1.50 

-1.09 ±0.02 

-1.31 ±0.18 

Holmberg II 

41.83 

41.72 

12.91±0.07 

12.80±0.03 

0.12 

0.07 

0 

41.86 

41.72 

-1.31 

-1.12 ± 0.02 

-1.12 ±0.17 

NGC 4214 

41.86 

41.82 

12.35 ±0.06 

12.06±0.03 

0.29 

0.78 

0.49 

42.17 

42.02 

-1.00 

-1.12 ±0.02 

-0.81 ±0.17 

NGC 5253 

41.87 

42.91 

12.87±0.07 

12.40±0.03 

0.48 

1.63 

1.14 

42.53 

42.36 

-0.64 

-1.13 ±0.02 

-0.46 ±0.17 

NGC 1569 

40.83 

40.94 

15.43 ±0.23 

14.76±0.09 

0.67 

>2.91 

>2.19 

41.99 

41.82 

-1.18 

-0.66 ±0.01 

-0.99 ±0.17 

NGC 4449 

42.70 

42.61 

11.20 ±0.03 

11.05±0.02 

0.15 

0.74 

0.53 

43.00 

42.82 

-0.17 

-0.41 ±0.01 

0.02 ±0.17 





Post-Starburst Galaxies 







Antlia 

39.02 

39.27 

17.89±0.71 

16.88±0.25 

1.01 



39.02 

39.27 

-4.15 

-3.96 ±0.15 

-3.97 ±0.33 

UGC 9128 

40.03 

40.03 

16.48±0.37 

16.08±0.17 

0.39 



40.03 

40.03 

-3.14 

-2.76 ±0.06 

-2.95 ±0.23 

NGC 4163 

40.61 

40.63 

15.69±0.26 

15.24±0.12 

0.44 

0 

0 

40.61 

40.64 

-2.55 

-2.40 ± 0.04 

-2.37 ±0.20 

NGC 6789 

40.29 

40.35 

16.89±0.45 

16.34±0.19 

0.55 

>0 

>0 

40.29 

40.35 

-2.88 

-2.39±0..02 

-2.69 ±0.25 

NGC 625 

41.43 

41.41 

14.22±0.13 

13.89±0.06 

0.34 

1.21 

0.85 

41.91 

41.75 

-1.26 

-1.56 ±0.02 

-1.07 ±0.17 


Note. — Column 1— Galaxy name. Columns 2 — 3 log of FUV and NUV luminosities calculated from the fluxes in Tableusing Equations and Columns 4 — 6 GALEX FUV and NUV magnitudes and FUV-NUV 
colors based on the measured fluxes with no corrections for extinction. Columns 7 — 8 FUV extinction based on TIR fluxes esti mate d from the measured Spitzer MIPS fluxes reported in[Mc^uinn_et_^IJ {201^ using the 
model of the TIR/FUV flux ratio from[Buat_et_a^ (200^) and the empirical relation given by IDale &: rfeT^ 120031 . See Sectionfor a full description. The UV luminosities for Antlia and UGC 9128 were not adjusted 
as both systems were not detected in any of the MIPS images. Column 9 — 10 GALEX FUV and NUV luminosities based on the measured measured IMcQuinn et alJ I201^ and corrected for extinction using the values 
in Columns 4 and 5. Column 11— SFRs calculated using the relation from [Hao_et_a^ (20^^) and the extinction corrected L^puv Column 9. Column 12— Average SFRs over the last 100 Myr derived from resolved 
stellar populations using CMD fitting techniques. Column 12— SFRs calculated using the relation from this study (Equation Ilol and the extinction corrected L^p^V Column 9. 
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Table 4 

Comparison of Measured UV Luminosities with Predictions 


Galaxy 

(1) 

GALEX 

T 0 

72T 

PEGASE/ 

GALEX 

(3) 

FUV 

MATCH/ SLUG/ 
GALEX GALEX 
(4) (5) 

GALAXEV/ 

GALEX 

(6) 

GALEX 

T 0 

777 

PEGASE/ 

GALEX 

(8) 

NUV 

MATCH/ SLUG/ 
GALEX GALEX 
(9) (10) 

GALAXEV/ 

GALEX 

(11) 





Starburst Galaxies 





UGC 4483 

3.25E-I-40 

0.69 

1.91 

1.31 

1.17 

2.41E+40 

1.28 

1.38 

1.19 

0.85 

UGC 6456 

1.02E+41 

0.37 

1.38 

1.23 

0.82 

8.01E+40 

0.66 

0.97 

1.05 

0.56 

DDO 165 

2.48E-I-41 

0.60 

1.99 

2.22 

3.41 

2.68E+41 

0.82 

0.98 

1.38 

1.68 

IC 4462 






5.43E+41 

0.58 

0.41 

0.70 

0.49 

NGC 6822 

7.13E+39 

1.97 

4.75 

2.69 

3.41 

4.19E+39 

5.04 

6.06 

3.09 

3.47 

NGC 4068 

2.83E+41 

0.43 

1.18 

1.45 

1.83 

2.12E+41 

0.81 

0.92 

1.30 

1.27 

NGC 2366 

8.40E+41 

0.64 

0.97 

0.94 

0.69 

5.81E+41 

1.23 

0.78 

0.92 

0.60 

ESO 154-023 

5.94E-I-41 

0.36 

1.01 

1.44 

1.93 

4.62E+41 

0.69 

0.70 

1.24 

1.38 

NGC 784 

4.66E+41 

0.45 

1.29 

1.93 

2.64 

3.86E+41 

0.84 

0.88 

1.57 

1.77 

Holmberg II 

7.27E-I-41 

0.46 

1.45 

1.26 

1.08 

5.31E+41 

0.88 

1.09 

1.16 

0.85 

NGC 42114 

1.48E+42 

0.31 

0.39 

1.05 

0.96 

1.04E+42 

0.61 

0.36 

1.00 

0.80 

NGC 5253 

3.35E-I-42 

0.32 

0.22 

0.49 

0.37 

2.31E+42 

0.63 

0.20 

0.48 

0.30 

NGC 1569 

9.80E-I-41 

0.47 

1.26 

1.54 

3.80 

6.57E+41 

1.03 

1.20 

1.55 

3.30 

NGC 4449 

l.OOE+43 

0.58 

0.37 

0.75 

0.64 

6.60E+42 

1.16 

0.34 

0.77 

0.57 





Post- 

•Starburst Galaxies 





Antlia 

1.05E-I-39 

1.21 

1.80 

1.65 

0.70 

1.86E+39 

0.97 

0.67 

0.63 

0.23 

UGC 9128 

1.07E+40 

0.46 

1.82 

1.87 

1.58 

1.08E+40 

0.79 

1.17 

1.26 

0.91 

NGC 4163 

4.12E-I-40 

0.26 

1.14 

1.12 

1.52 

4.33E+40 

0.42 

0.65 

0.72 

0.75 

NGC 6789 

1.96E-I-40 

0.49 

1.71 

2.10 

4.65 

2.26E+40 

0.74 

0.85 

1.23 

2.03 

NGC 625 

8.19E+41 

0.07 

0.16 

0.34 

0.72 

5.57E+41 

0.17 

0.15 

0.34 

0.58 


Note. — Column 1 —Galaxy name. Column 2—GALEX FUV band luminosity (ergs s based on the FUV fluxes measured in the equivalent HST field 
of view and corrected for extinction. Column 3—6 Ratio of the predicted FUV luminosity to the measured FUV luminosity based on SFRs(t,Z) derived 
from optically resolve d stellar popula tion and the synthetic spectra l model PEGASE 2.0 dFioc Sz Rocca-Volmerangelfl997n . synthetic stellar populations 
generated from MAT CH (|Dplphi^l2002D . the stochastic model SLUG ((d a Silva et al ] [20Ta . and the synthetic galaxy spectrum from the GALAXEV code 
jBruzual & Charlotll200^ . Columns 7—11 are the same as Columns 2 — 6 but for the GALEX NUV band pass. Note the luminosities and predictions 
listed for NGC 6822 are the average luminosities over three fields of view. 


Table 5 

Stellar Atmosphere Models 


Model 


Stellar Library References 


Starburst99 

Pegase 

MATCH 

SLUG 

GALAXEV 


^ejeun^^^^^ ^^97,^98) for plan-parallel atmospheres which includes models from jGirucz ^^9^ 

Hillier Mi ler (19981 atmospheres for stars with strong winds 
£auldracl^^al^ (20^) models for O star atmospheres 

IGiru^ ^^9^ stellar spectra compiled by ^ejeun^^^al^ ^^97,^9^ for plane-parallel atmospheres 
^asteU^^MGirucz (200^ 

Black-body spectrum for non-white dwarf stars with T> 50000 K 

Leieun^^^aL ^^97,^98) for plane-parallel atmospheres based on the original model spectra fromJGirura (199^) 
^^^^^^^^^^^^^implementation of atmospheres for O stars from JLllie^^^^ll^ ^^9^ 

Maeder (1990) for stellar winds 

^ejeun^e^al^ ^^97,^98) for plane-parallel atmospheres which includes models from j<uru^ ^^92) 

Weste^ I -'III 11 r ^ester^^^a^ ^2002) for semi-empirical calibrations of model-atmospheres at non-solar metallicities 































